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ABSTRACT 


The anomalous microwave radiation from interstellar OH has been the 
subject of intensive investigation since its discovery in 1963. An 
investigation of the mechanisms put forward to explain the microwave 
radiation now indicates that several different far- infrared .lines with 
.wavelengths between 50 and 200 microns can also show the anomalous maser 
emission. Further, different pumping mechanisms produce different far- 
infrared lines, so that spectral studies in the far-infrared to be made 
in the near future will provide additional information on the conditions 
and pumping mechanisms of the maser regions. 

The energy structure and transition rates of the OH radical have 
- - ' 

.been carefully studied in rotational and vibrational energies to provide 
the basis for studies of the population distributions under various ex- 
citing mechanisms. The absorption of ultra-violet or near infrared radia- 
tion which Litvak and others have proposed to explain the microwave 
emission has been studied, as well as the absorption of low temperature 
blackbody radiation and collisions with molecular hydrogen. Each of 
these provides population inversion in some infrared lines, while dipole 
collisions with charged particles does not. 

In applying the results of the population studies to the transfer 
equation a model of a uniform spherical cloud of OH is used with the 
pumping process occurring uniformly throughout the region . For high 
pumping rates the emission is unsaturated and the intensity is an exponen- 
tial function of the size of the cloud-. For lower pumping rates or 
larger sizes the region of unsaturated emission shrinks until the entire 
region is saturated, the intensity then becoming a linear function of the 
size of the emitting cloud. 
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— The population solutions and the intensity formula are combined to 
show that under the various exciting mechanisms covered far infrared 
line radiation between 50 and 200 microns in wavelength will be emitted 
isotrppically from the OH clouds at intensities detectable by a 36 inch 
airborne infrared telescope to become operational next year at the NASA 
Ames Research Center. The solutions of the transfer equation rniay also 
be applied to the microwave maser emission observed from interstellar 
water vapor and other organic molecules to provide information on the 
conditions in the many dense molecular clouds which have been recently 
discovered. The line widths of the emission and the observed sizes of 
the emitting clouds is briefly considered. 
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CHAPTER I 


INTRODUCTION 

Interstellar gases are largely composed of atoms ,of ions at a density 

of about one atom per cubic centimeter. The ultraviolet radiation from 

stars can ionize such gases, producing electrons of temperatures charac- 

4 o 

teristic of the stellar radiation, about 10 K. There are also sub- 
stantial amounts of neutral hydrogen and even moleculai* hydrogen. The 
gas apparently has two stable states: The first is a dilute gas of about 
0.1 atoms per cubic centimeter with a high temperature and a high degree 
of ionization; the other state is a much denser, cooler gas cloud with 
large amounts of dust particles. 

When a gas cloud becomes dense enough collisions can provide an 
important cooling mechanism. This will overcome the heating by ultra- 
violet radiation, allowing the gas to contract and cool further. This 
in turn increases the dust concentration and the opacity, and with the 
shrinking size of the cloud less and less heating radiation is absorbed. 
The dense cloud of gas and dust apparently reaches a low temperature 

4 6 

stable state with densities near 10 to 10 atoms per cubic centimeter, 
and temperatures of ten to a few hundred degrees K. 

These dense dust clouds have only recently received attention, 

partly because their darkness makes most optical observing techniques • 

* 

useless. Microwave and infrared radiation can easily penetrate these 
clouds because of the small size of the dust particles and recent advances 
in microwave and infrared techniques have revealed some startling charac- 
teristics in dense clouds. They represent a significant fraction of the • 
mass of our galaxy and are important in many galactic phenomena including 
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the spiral structure and star formation. Most startlingly, it is in 
these dark clouds that many complex molecules have been discovered. 
Molecules were first found in interstellar space in the 19 ^ 0 ' s 

when they were identified in optical spectra taken by well established 

* *i* 

techniques. Three molecular free radicals were found: CH, CH J , , and CN. 

The development of radio astronomy led to the discovery of the next mole- 


cule in 1963? the OH radical (Weinneb et al . 1963) . These sparse dis- 

coveries and the concept that the interstellar gas is of low density and 
high ionization led to the expectation that molecules in interstellar space 
would be both rare and simple. It seemed reasonable however that dark clouds 
would be the most likely place for molecules because of their higher density 
and lower temperatures .- Furthermore, the dust should provide a shield 
against the ionizing radiation of background stars and in 1968 a search of 
such clouds revealed the first polyatomic molecules, NH^ and H^O (Cheung 
et al. 1968). Both of these have convenient microwave spectra with wave- 
lengths long enough to penetrate dense clouds easily. The organic molecule 
HgCO was discovered soon afterwards, and since then the number of molecules 
found in interstellar space has increased rapidly. The total number of 
molecules and isotopes discovered has reached nearly thirty, and it seems 
that many others will be found as further searches are made with' increasing 
sensitivity in microwave and infrared spectral regions. 

The microwave lines of the OH radical have been observed in the 
absorption and emission spectra of several hundred galactic regions. 

The OH seen in absorption appears to be associated with regions of neutral 
hydrogen as suggested by line widths and radial velocities which agree 
closely with the observations of 21 -cm hydrogen lines . In a very few 
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cases the four ground state OH lines are observed with intensities charac- 
teristic of local thermodynamic equilibrium, and comparison with the 

_ ii 

hydrogen observations indicates an abundance ratio of OH/H between 10 
-7 

and 10 , the larger value being associated with the galactic center 

region . 

For most of the sources, neither "the OH or the related H^O have energy 
levels populated according to the Boltzmann distribution. This was 
indicated for OH shortly after its discovery in absorption by line 
strength ratios which differed significantly from the LTE predictions. 

The degree of anomalous excitation indicated by the observations of OH 
in emission was even more unexpected. The OH emission is now believed 
to result from maser amplification, but in the eight years since its 
discovery the most notable properties of the OH emission remain unexplained; 
even the mechanism which populates the energy levels is not yet known 
with any certainty. Several mechanisms have been proposed for populating 
the energy levels and this paper considers the resulting infrared line 
emission for the most interesting of these as a potential observational test 
for the relevant mechanisms. 

A. Observations of OH Emission 

Anomalous OH emission was discovered close to HI I regions in 19 65 
by Gundermann (1965) and by Weaver et al. (1965). This emission was • 
strongest at 1665 MHz and was highly polarized, usually circularly. These 
properties along with very small angular sizes as established by inter- 
ferometry applied to all seven sources found in the following year or so. 
These early searches were usually made only in the two main lines and 
were confined to the brightest HII regions. In 1966-67, Goss (1968) 
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extended survey work to weaker HII regions and to galactic non-thermal 
continuum sources, as well as searching in all four ground state lines. 

More sources of the main line type were found as well as several sources 

of a new type which emit only at 1720 MHz and which are not associated 

/ 

.with HII regions. This latter type of emission was subsequently found 

k « 

to be unpolarized. The third major class. of. emission sources which radiate 
predominantly at l6l2 MHz was discovered by Wilson et al. ( 1970 ) to be 
associated with stellar like objects whose continuum radiation is primarily 
infrared. Finally, the first apparently normal OH emission of very weak 
intensity was found in dust clouds by Heiles ( 1968 ). By 1 969 , thirty 
anomalous OH emission sources had been discovered from the northern hemi- 
sphere and nine more from the southern hemisphere, all of which have 
been well studied. More recently a large number of new OH emission sources 
have been found associated with continuum sources (Downes 1970; Turner 
1970a). Turner (1970b) proposed the following classification scheme. 

B . Classification of OH Emission • 

Class I. These sources emit very strongly in one or both of the 
main lines, with the 1665 MHz line nearly always the stronger, and weakly 
or not at all in the satellite lines. All four lines in this class of 
emission sources are highly circularly polarized, with the degree of 
circular polarization usually being stronger for the most intense sources. 
Some of these sources have features that are linearly polarized to a degree 
higher than can be explained by resonant scattering, with no circular com- 
ponent. The emission profiles show little similarity among the four OH lines. 

* / 

They consist typically of many narrow components of widths less than 2 kHz. 
The line ratios in these sources vary strongly with velocity and time and 
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do not resemble LTE predictions; the time scales are between nine days and 

two months. If due to thermal broadening the linewidths correspond to 

o -O 

temperatures below 50 K with the narrowest giving a temperature of o K, 

Very long Baseline Interferometry indicates that Class I sources with 

several radial velocity features are close collections of separate sources, 

with linear sizes about 10 AU for- individual features -and separations of 

about 100 AU for features of a single source. Continuum maps of these 

regions indicate that the OH emission is near the continuum centers but 

does not coincide. There is general evidence to associate Class I emission" 

/ 

with infrared sources (Eliasson and Bartlett 1969 ; Wilson et al. 1970), 
and a few are associated with compact HII regions (Mezger et al. 1967 ) • 

The most -heavily studied of the Class I sources are in W3, Orion A, W 5 I 
a and b, DR21, Sgr B2, W^9, and NGC 633 k. 

Class Ila. These sources show anomalous emission only in the 
1720 MHz satellite line, anomalously strong absorption in the l6l2 MHz 
line at the same radial velocity, and normal absorption in the main lines. 
These sources are not 'highly polarized, with typical values of 10$ circu- 
lar or linear polarization in the emission line and none in the absorp- 
tion lines. These sources are speculatively associated with the shells 
of supernova remnants, although the distances are not known. No high 
resolution continuum mapping has yet been \ione on these because of the' 
weakness of the emission. Typical members of this class are W28N, wMt, 

CasA and NGC 6857 . 

Class lib. This group of objects emits most strongly in the l6l2 MHz 
line with no emission and little if any absorption .in the 1720 MHz line. 

The two main lines show emission much weaker than either the l6l2 MHz line 
or the main lines of Class I sources. The intensities in the l6l2 MHZ 
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line is much, stronger than the Class Ila emission and occasionally stronger 
than the brightest Class X objects. The l6l2 MHz linewidths are typically 
15 kHz with 1665 MHz linewidths about 3 kHz in the same feature. The 
strong satellite line is nearly unpolarized in all objects of this class, 
and the weaker main lines are circularly polarized to as much as 65$. All 
of the sources in this class are associated with IR sources listed in 
the CIT catalog to the resolution of VLBI measurements, with several small 
spots of individual sizes around O'.'l in a region about the IR star of 2" 
extent. Objects in this class include VYCa Maj, NML Qyg, NML Tau, CIT 3 > 
and CIT 7- 

Class III. Heiles (1968) discovered weak OH emission with LTE line 

ratios in a number of dark dust clouds. Both the clouds and the OH 

emission are extended,, up to 2° in one case. There is no continuum radio 

emission and little or no 21 cm hydrogen emission. The excitation tempera- 
o ■ +6 

tures are about 5 K with densities greater than 10 molecules per cubic 
centimeter. There is no apparent polarization in any of the lines in 
these sources, typical of the thermal characteristics. Careful observa- 
tions by Turner and Heiles (1970) indicate some anomalies in the 1612/1720 
satellite line ratio similar to, but weaker than, the Class Ila objects. 
There is some speculation (Booth 1972) that this class is a weak case of 
Class Ila or b. 

C. Characteristics of Two Strong OH Sources 

The largest infrared telescopes available for observations beyond 
100 microns are less than 1 meter in aperture, so that sensitivity -is 
rather poor in comparison with radio det.ectors. The far infrared radiation » 
detectable from the OH molecule will probably come from the strongest 
OH radio sources. Any reasonable pumping mechanism for the ground state 
maser emission will involve infrared transitions so that the infrared 
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brightness temperatures will be related to the microwave emission by the 

ratio of the wavelengths, hence about 10 8 K. The limited sensitivity of 

the IR telescopes restricts detectability to the brightest of the OH 

sources. For this reason, we will consider modeling the two brightest 

class I sources, Orion A and W3, in later chapters. 

The W3 source consists of several spots of emission about 8 AU in size 

in a circular arc about 3000 ATJ across. Microwave emission has been 

2 

detected in the ground state and in four excited states, ' ^ 3 / 2 ^ ~ 

3 / 2 , 5 / 2 , 7 / 2 , and = 1/2 and 5/2. The relative intensities in 

these lines are 1.00, 0.03, 0.06, and 0.01, 0.0005 respectively. The 
' 2 

emission in the ~ 5/2 lines may be a false detection (Turner 1972) 

W3 is associated with a compact HII region as well as a point - inf rare'd 

source. The brightness ^temperature of the 1665 MHz line is about 10 J K 

and the linewidth about 0.8 kHz. There is also an anomalous H^O emission 

source nearby. Time variations are seen of the order of a few weeks in 

the 1665 and 1667 MHz lines, but with no correlation between them. The 

polarization remains constant with time during these intensity variations. 

28 

The total intensity in the lines assuming spherical emission is about 10 
ergs per second-. 


The Orion A source is apparently a single small spot about 3 X 10 


,13 


centimeters across with a nearby IR point source and a large background 

12 o 

IR nebula. The brightness temperature is about 10 K and the linewidth 

10 kHz. There is an H^O emission source correlated with the OH source 

• 1^ o 

which has a brightness- temperature greater than 10 K. 
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D. Summary 


In the following chapters we will discuss in detail the energy 

l6 1 

structure of the OH molecule, specifically the 0 H isotope, and cal- 
culate the transition rates.- This information will be applied to the 
pumping mechanisms to predict the population distributions for varying 
parameters of the OH regions. The equation of transfer for a simplified 
spherical region will be solved for various degrees of saturation and 
analyzed for relevant emission characteristics. Finally, the two results 
will be combined to predict the expected detectable infrared emission for 
some typical OH sources and to describe the relevant physical conditions 
in the OH regions. 
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CHAPTER II 


THE MATHEMATICAL FORMULATION FOR OH MASERS 


A . Notation of the Energy States 

The energy structure of the OH molecule is as complicated as any 

* 

diatomic molecule, with contributions from the electronic excitation, 
the internuclear vibration, the molecular rotation, the electron spin, 
the nuclear spin, and the interaction of each of these with the others. 
There are no less than nineteen quantum numbers appropriate to the 
description of the energy levels of the two electronic states of 
interest in the pumping mechanisms of OH masers. The following descrip- 
tion of these quantum numbers is a useful adjunct to the understanding 
of the energy structure and the emission which accompanies an OH maser. 

The coupling between the two nuclei of the molecule can be described 
as a nonlinear spring, or quantum mechanically as a particular potential 
function of the separation. Vibration of the molecule in this potential 
is quantized, and described by the vibration number V. This is the 
only energy in the molecule which has no atomic equivalent. 

The motion of the electrons in an atom takes place in a spherically 
symmetric field, hence the electronic orbital angular momentum L is a 
constant of the motion. In a diatomic molecule the symmetry of the field 
is only axial and as a consequence only this component of the orbital 


angular momentum is a constant of the motion. A precession of L occurs 
about the internuclear axis with a component along the axis of which 

takes the allowed values of L, L-l, L~2, . . .', -L. In the absence of 


perturbations the states with 




values of opposite sign, have the same 


energy, thus the component of the electronic orbital angular momentum is: 


A ■= l»J • 
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Note that except for A = 0 the states are twofold degenerate. In a 
manner similar to atomic states, these molecular states are designated 
as 2, n, A, <j>, ... according to A = 0, 1, 2, 3? • * • respectively. 

The wavefunction of the molecule in these axially degenerate states 
may be written in a manner which expresses the symmetry characteristics 
of the molecule: 



Y 


e 


= x e iAtp - x e - 1A<p 


- x e iA<p ♦ x e - 1A<f 


(II-A-1) 


" “ 

hence the degenerate states may be written as 2 , £ , II , n > '• • • • The 
separation of the symmetry states by A-doubling perturbations will be 
discussed later. 

A closer investigation of the structure of electronic states shows 
a multiple! structure due to electron spin just as in atomic states. 

The spins of the individual electrons form a resultant S with the 
quantum number S being half-integral for molecules with an odd number 
of electrons such as OH. For A = 0 the resultant S is fixed in space 
as long as the molecule is not rotating. In A ^ 0 states the orbital 
motion of the electrons gives rise to an internal magnetic field which 

i 

causes a precession of S about the internuclear axis with a component 
2 along the axis with the values S, S-l, S-2, ..., -S. 2 is not a 

well defined quantum number for A = 0. 

The total electronic angular momentum about the internuclear axis 
is the combination of orbital motion and spin. The quantized component 
of that total along the axis is 0 = |A + 2| . For A ^ 0 the electronic 
state splits into 2S +1 levels. Additional angular momentum in the 
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molecule is contained in gross rotation of the molecule and in nuclear 
spin. The nuclear spin quantizes just as the electron spin, -with resul- 
tant JE and component . The nuclear spin splits the states into 
21+1 further substates. The rotation of the molecule occurs perpen- 
dicular to the axis, with a resultant 0. Table 1 summarizes the quantum 
numbers in diatomic molecules. 

The several forms of angular momenta in the molecule' can interact 
in different ways, Hund first studied the various coupling schemes in 
detail, and his cases (a) and (b) are of interest (Herzberg I950) . In 
the OH molecule, the ground state (A = 1, a II state) is intermediate 
between case (a) and case (b), and the first excited state (A = 0, a S 
state) is a pure case (b) . 

In Hund's case. (a) 'the interaction of the molecular rotation with 
the electronic motion is weak, and the electronic motion is coupled 
strongly to the internuclear axis. The vector diagram, for this is: 



- \ 1 


The quantum number for the rotation can take -the values Q, Qfl, Qt-2j ... . 
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TABLE 1. THE QUANTUM NUMBERS FOR DIATOMIC MOLECULES 


Nuclear spin 

projection on molecular axis 

projection on space axis 

Electron spin 

projection on molecular axis 

projection on space axis 

Electron orbital angular momentum 

projection on molecular axis 

projection on space axis 

Molecular rotation 

Total orbital angular momentum . . . • 

projection on molecular axis . 

Total angular momentum without nuclear spin. . 
projection on molecular axis (absolute value) 

project ioii on space axis . 

Total angular momentum 

projection on molecular axis 

\ 

X 

projection on space axis .... 1 


I 



A 

\ 

0 

N 

K 

J 



Vibration 


V 


P P 



In Hund's case (b) the spin is only weakly coupled to the inter- 
nuclear axis, hence the electronic orbital angular momentum and the 
rotation form a resultant N and its quantum number can take the values 
A, A+l, A+2, ... . The angular momentum N and the spin S then coupl 
to form the total" angular momentum J with quantized values Nl-S, 

NtS-1, ..., N-S . The vector, diagram for this case is: 



B . Calculation of the Energy Levels 

As mentioned above the energy of an OH molecule is in several forms 
Since it is known that the energies associated with the interactions of 
the various motions of the molecule are much smaller than the energies 
of those motions, those interactions can be considered as perturbations 
in deriving the expressions for the energy structure of the OH molecule. 

In all cases the energy expressions in this paper will be given in 
term values of cm ^ so that the energy differences will be directly in 
units of inverse wavelength. To determine the energy in more common 
units the expressions should .be multiplied by Planck's constant and by 
the speed of light (h x c ) . 
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The energy states of interest in the pumping mechanisms are the rota- 
tion levels of the ground vibration state in the first excited electronic 
state, that is the J levels of .the V = 0 state of A = 0 (the 2 
electronic state) ; and the vibration and rotation levels of the ground 
electronic state, that is the J and V levels of A = 1 (the n state) . 
The 2 doublets of opposite symmetry will be considered, as will the jj 
levels split by symmetry, electron spin, and nuclear spin. Figure 1- 
diagrams the labeling and order of the states for A = 0, and Figure 2 
shows the typical pattern of the A = 1 ground state for each vibration 

level. The levels in the scheme are further split by the action of 
• \ 

the unpaired nuclear spin which is shown in Figure 3 for a typical J 
doublet . 

The vibration energy of a diatomic molecule is easily separated from 
the other forms of energy. With the vibration quantized by V = 0, 1, 

2, ... the vibration term value is: 



(II-B-1) 


where the constants oo e , x t y Q are determined experimentally for each 
electronic state, and higher order terms are ignored. The vibration 
may interact with the electronic motion of the molecule by altering the 
average moment of inertia of the molecule. This is taken into account 
by writing the rotational constants of the molecule as functions of the 
vibration number; 


-B^ 1 ) = B + a 
e e e 


= D e + Pe( V + l) ■ 


( v - 1) - + l) ! 

4b3 


°e ~ 2 

CD 


i4 


(II-B-2) 



Figure 1. The 2 energy levels 


K * J 

O 1 7/2 - 

3 ' 2 5/2 - 

P 1 5/2 + 

2 3/2 +* 

, r 3/2 - 

2 1/2 - 

0 1 1/2 + 


Figure 2. The II energy levels for each vibration state 


K * J 

4 1 9 /2 ^ 

3 £ 7/2 _ 


2 

1 


1 ' 

1 

1 ' 

1 


5/2 + 

4 - 

3/2 


n 3/2 states 

Figure 3* The hyperfine structure 




states 


F= 

2J + 1 

2J - 1 

2J + 1 
2J ~ 1 


* J 

2 / 7/2 

5/2 

3/2 


*This is the label of the spin splitting. The prime refers to the 
A-doublets . 


15 



where higher order terms are again ignored. The subscripts used in 
( II-B-l) and (lI-B-2) refer to the electronic state being considered. 

The letter p will be used for the n states and s for the £ states 
throughout this paper. 

The rotational and electronic energies of the molecule take a more 
complicated form. After Van Vleck (1929) the hamiltonian of the molecule 
including the electronic and rotational terms and the L • S interaction 
may be written: 


H = B[J(J+1)-A 2 ] + BS(S+1) + Ah z S z - 2BJ • S + B(L x 2 +L y 2 ) 


+ (2BfA)(L x S x +L y S y ) - 2B(J x L x +J y L y ) . (II-B-3) 


Here the z direction is along the internuclear axis, A and B are the 

fine structure interaction and rotation constants, and the quantum numbers 

2 2 

have been described. The term B(L x +L y ) is nearly constant and will 
be ignored since it only shifts the energy by a constant amount (see also 
Fontana 1962; Hill and Van Vleck 1928) . 

The eigenvalues of the wave equation with the above hamiltonian are 
needed. There are six states to be considered: IIy 2 > U_y 2 ’ II 1 / 2 > R-1/2’ 

£ 1/2 3 E 1/2’ T ^ e iHteractions between all pairs of states form a six-by- 
six determinent that factors into two three-by-three determinents by 
using the symmetric-antisymmetric wavefunctions mentioned previously. 

Using the matrix elements of Van Vlteck one determinent is: 


k-X p I] 

* 

p 0 -X £ 

H & . y-\ 


(iI-B-ii) 


% 
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■with: 


k = a+ 6 , [1=0+5 (II-B-5) 

and the matrix elements are listed in Table 2. The other cubic determinent 
is identical except: 

k = a - 6 , [ 1 = 6 - 5 . (II-B- 6 ) 

The eigenvalues may be found by diagnolizing the above to successively 
higher order (Dousmanis et al . 1955 )* 

The zero order terms are the energies of the state with no 
interactions ([i = = 0): 



1/2 (II-B-7) 

4^3 = ± l + 4&2 ] + | • 

The zero order energies then are:’ 

= E + B [J +• 1/2] [J + 1/2 ±'l] 

s ± 

1 ' • (II-B-8) 

E^ 0) = B (J + V2)(J -.1/2) ±% B \jk(j + 1/2) 2 + \(X~b) : 

^ 1 / 2 , 3/2 P ^ P 

Note that to this order, not even the symmetry differences in the T[ state 
are split. The interaction between the 2 and jt states in first and 
second order is needed to describe this splitting. 

We proceed by expanding the determinent ( 5 ) and writing the eigen- 
values as perturbations of the zero order terms above: 
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TABLE 2. THE MATRIX ELEMENTS OF THE HAMILTONIAN 
(Van Vleck 1929 ) 

a = B [J(J+1) + 1/k ] + E 

s 

P = B p [j(j-1) + 1/4] - | A 

Y = B p [J(J+l) - 1/k] + | A' 

6 = B (J + 1/2) 
s 

ft = B p [(J - 1/2) (J + y2)] l/2 
5 = 2(j + l/2)<n|BL y |s) 

H = 2[(J - 1/2) (J + 3/2)] 1 / 2 (nK y |s) 
e= (n|AL y + 2 HL y |s) 


18 



, * * * * * * 
k£v + fjiRTi + n n a - pim - ypu - k ee 


- (kp + ky + PY - w* - UP - Till )X + (k + p + y)X 2 -X 3 


= (x - XjKx - x 2 )(x - x 3 ) 


(II-B-9) 


where: 


\ = 4°) + x[ X) , i = l,2 j3 . 


(II-B-10) 


If the \ from (II-B-10) are inserted in (II-B- 9 ) and terms higher than 
first order in X^^ are omitted, these first order corrections to 
(II-B-7) are found: 


(1) , * *. , 

’ = (wx + Tin )/E 




( 1 ) (mm. +Tffl ) T 1 

^2,3 ' 2E 


EB X 
P 


(P-Y)-(W* -HT*) + + ^*1)] 


(II-B-II) 


where 


X = +[4(J + 1/2) 2 + (y-4)y] V2 y = A/B p . 


(II-B-12) 


This procedure is repeated, writing: 



4 0) 


♦ *< 2 > 


(lI-B-13) 


for the second order corrections. Using this in equation (ii-B~9) gives 
the complex expressions for the perturbations: 
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( 2 ) 


- hi * 


v - 0 ,114 p + mn y + "Hu 

E 


* * / * *\ / \" 
a - nm + nu JloHEts) 


. ( 2) = ± 1 

2 ’3 e 2 b X 


/ * *- 

HU + THl 


)M-( 


* * * 

HU P + T1T1 Y + U&'H + tlH & 


'P*Y ± B X 


P \ / ,, \ / ^ ^ \ 

1 - (ai-E+6)[p'p + YHi-i “ M.&TI - R & TU 




* *x/^ ±B p x 


T1T| + HU 


B 


1 r / \2 * * * / 2 *2 2 *oV 

^2 [(P-y) MU THl + £& [H U + Tfr 2 j 


* * * / \ / * *\/ 
- 2(iji Till CS - (P-Y HU ~ TIT] Hi 


* *\/ * * *\ 
HH - Till J( uSTl + Tlu 6 I 


( II-B-14) 


The upper of the double- signs in equation ( II-B-14) corresponds to the 
Il-jyg state and the lower to the n^y 2 sta-fce of the OH molecule. 

Evaluating the energy terms of equations (II-B~7)> (H~B-ll), and 
( II-B-14) with the matrix elements of Table 2 gives the energy levels of 
the molecule with giving the energy of the E states and and 

X. giving the energy of the ^/2 anc * ^3/2 sta ^ es • For the H states, 
the energy is written: 


w = Wi ± - w 2 


(II-B-15) 


where \V^ explicitly is the A-doubling energy splitting, may be 

expressed as a sum of first and second order terms: 


w = wjj 1 ) + 

2 A A 


(II-B-16) 
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The expressions for the energy of the j] states are: 


B^X 


\ = B^(J - 1/2) (J + 3/2) ± -|- + ^ 



B 


B 


1 + ? (j-1/2) (j+3/2) - ^ (J4-1/2) 2 - (j + l/2) 


B X 
2\ . P 


2 . 


E 




4b 

s 

E 


(j+1/2) 2 (j-1/2) (j+3/2 ] 


W + i 2 ( J - 1 ''2)(j + 3/2) - ^ (^f 


_ fJ-1/2) (J+3/2) 
B 

P 


a. 


m 


- 2a B (2-\)[J(J+l)+3A] + 2 

p p i6p 


j 2 + J(j+i)+y4j + p p 2 (j+i/2) 2 [( 2 -x) 2 +i] 

« p 3 (2-X)-j 


P t> ( V B s ) 2 ■ 2 • 

-P--g— S - J 2 (J+1) 2 


(XI-B-17) 


where 


B' = B -p (l - ^ 
P P P\ Ei 


a. 


a 1 + 
p\ 


B -B \ 

p s\ 


5p;)( 2E 


C„ = 


V (1 


B s B 1 - 
P P\ 


XB B -B 
P JP S 

2E E 


2 2 

-SpTT- 


(lI-B-18) 


\ 11 

( , XB p 

B -B \ 
P s 1 

) 2 P p 1 

^ + 2E 

" -*5e-| 


a b 
p s 

2E 


and 
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W A = _a p (j + 1/£ )( j:X ±X X ] ( J _1 /2) ( J+l/2) ( J+3/2) (II-B-I9) 


and 


p (J+1/2)ol 

w ± A 

A ex 


2 + B 

P\ 2 


2 (j-l/2)(j+3/2) ± x)j (B s -B p ) (j+l/ 2 ) 

- 2B s ^(j+l/2) 2 ^j(j+l)-7A+x/2j - (j-1/2) (j+3/2) j 


(j-1/2 



, ( J+1 /2) n 
EX 


^2 (j-1/2) (j+3/2) ±,xj^ 2 B g (j(j+l)-lA) - 2 B p ( J-1/2) (j+ 3 / 2 ) j 


B s ,(2(j + 1/ 2) 2 -x) 


( j+l/2) a p 
EXAp 


B S (X - 2 T X) 


C (J-1/2) ( J+l/2) ( J+3/2) 


B X~ 
P 


(II-B-20) 


where 


c 0 = -20 b (\ 2 -4\+6) + (30 2 + 4b 2 )(2-x) - a„ 3 /B . (n-B- 21 ) 


p- p \‘ 


p "p 


In equations (lI-B-17) to (II-B-21) the upper sign corresponds to the II ^ 
states of the OH molecule and the lower sign to the s ^ ates (Dousmanis 1955)* 

Equations (lI-B-7), (II-B-12), and (II-B-14) may be evaluated for the 
energies of the 2 state. For each rotational level with rotational 
quantum number N there are two states that correspond to the two 
orientations of £ with respect to N in the OH molecule. For these 
paired levels the energy is given by: 
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w = c + b'n(w-i) 

z ± 3 s 



a 2 b “* 

s p 



(N4-1/2) 




(V 2 ^ 



N(W-l)(Nfl/2) 


where: 




and 


B' 

s 


B s +S *. 



CL. 


B -B 
_P s 

2E 


a 


B +B 




E 


There is only one other effect that must be accounted for to model 
the energy of the OH molecule to the accuracy of the second order terms 
given above. The rotation of the molecule alters the effective equili- 
brium position of the mass-spring-mass model. This increases the 

\ 

moment of inertia with increasing rotation and can be expressed by modifying 
the rotation constant in addition to the expression (II-B-2) which corrects 
for the average moment of inertia due to vibration. We write: 


B< 2 > = B* 1 ) 
e e 



(II-B-23) 


The higher order terms in (II-B-3) may also be similarly corrected: 
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(II-B-24) 


a 


all 



P 


e 




f 

e 


For II states: 


f 

P 


= J 2 + 


t* X + 2(J+l/2)-TrT7p\ 

l Tx )( J+1/2 ) 


(II-B-25) 


For 2 states: 

f = N(Nfl) . (lI-B-26) 

s 

The upper signs in ( II-B-24) are for the Iljyg states and the lower ones 
for the n 3/2 states. • 

The expressions (HrB-22) give the energies of the 2 states to 

sufficient accuracy for our purpose. Expressions (II-B- 17 ) to (II-B-20) 

give the energy of the ]J states to the resolution of the A~doubling. 

The states are further split as mentioned by the action of the nuclear 

spin. To the accuracy of the above expressions, this hyperfine splitting 
* 

is given by W = W + where ■ ' 


\f = ±5x 


2a(± 2X + 2-X) + b 


4'(j-i/2)(j+3/2) ± X + 4-2xj 


+ 0(± X + h-2 X ) ± | fiij ( J+l/2) | 


(ii-b-27) 


where 
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( II-B-28) 



The upper sign in front of X in (lI-B-27) applies to the state 

“ » 

while the lower one applies to the state . Only the d term above 

will give rise to any splitting of the energy levels, since the upper 
sign here corresponds to the hyperfine level with F = J + 1/2 while 
the lower sign is for the F = J - 1/2 level'. The other terms only 
make the splitting asymmetrical about the A-doublet level. For transi- 
tions between the rotational levels the b term will shift the hyperfine 
lines slightly from the unsplit frequency. 

No further splitting of the II levels will occur without the 
application of an external field. A magnetic field will split each 
hyperfine level into 2F + 1 magnetic sublevels which can give rise to 
polarization, in the emission, but -these splittings and the hyperfine 
splittings are smaller than the. typical width of the far infrared lines, 
and certainly smaller than the resolving ability of current infrared specto- 


scopes . 


C. The Molecular Constants and the Energy Levels 

The spectrum of the OH molecule and its isotopes has been observed 

/ 

by several groups. A reasonably accurate fit of their data to the 
previous formulae for the energy levels gives the following set of 
molecular constants (Herzberg 195 O; Dousmanis et al. 1955 ; Dicke and 
Croswhite 1962): 
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Vh 

n 

S 


m e 

3735.21 cm -1 

3184.28 cm -1 

‘Ve 

82.81 " 

97.87 

it 

cu y 
e e 

— 

— 


B 

e 

18.871 ” ‘ 

(17.355) 

»» 

«e 

-0.7703 ” 

-0.807 

11 

^e 

0.0192 ” 

-O.OO83 

11 

E 

e 

0 

35.875 

It 

Y 

e 

■ . , * 
- 7.444 

— 


C 

e 

0.60 

— 


F 

e . 

1.20 

— 



The hyperfine constants that have been measured are: 

a = 0.002 cm 1 , 

d = 0.00190 cm -1 . 

-1 

The energy levels in cm resulting from this set of values and 
the previous formulae are listed in Table 3, where f^ and f^ are the 
A-doublet levels in the 11^ ^ state and f^ and f^ are the doublet 
levels in the Iljy^ state. K is given by K = J - 1/2 for the 
state and by K = J 4- 1/2 for the Iljy state. Table 4 gives the 

' The value of the intermediate coupling strength is not constant for 
all vibration levels. Rather it has been shown to vary, and the most 
commonly reported values are (Benedict et al . 1953): 
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TABLE 3. ROTATIONAL "AND VIBRATIONAL LEVELS OF THE II 
STATE OF OH ( IN cm” 1 } 


K 

f i 

5 

f 2 

f 2 

V = 0 

1 

00.00 

00.05 

126.43 

126.12 

2 

83.70 

83.90 

187.71 

187.47 

3 

201 .go ' 

202.37 

289.01 

288.83 

4 

355.09 

355.87 

429 • 45 

429.23 

5 

543.5^ 

544.82 

608.15 

608.16 

v - 1 

1 

3568.31 

3568.49 

3695.34 

3695.03 

2 

3649.10 

3649 .28 

3754.15 

3753.95 ■ 

3 

3763.03 

3763.40 

3851.38 

3851.10 

4 

3910.55 

3911-26 

3986.14 

3985.99 

5 

4091 .95 

4093.10 

4157.80 

4157.75 

v = 2 

1 

6971.02 

6971.I8 

7097.38 

7097.30 

2 

7048 .85 

7049.16 

7155.11 

7155.20 

3 

7158.60 

7159.00 

7248.20 

7248.02 

4 

7300.50 

7301 .24 

7377-48 

7377.32 

5 

7474.98 

7476.04 

. 7542.10 

7542.10 

.V = 3 

1 

10210 ..31 

10210.77 

10335.27 

10337.60 

2 

10284.31 

10284.40 

10390.21 

10392.22 

3 

10390.59 

10390.94 

10480.59 

10481.33 

4 

10526.99 

10527.5^ 

10605.22 

10605.07 

5 

10694.61 

10695.47 

10763.OO 

10763.07 


TABLE 4. ROTATIONAL LEVELS OF THE 2 STATE OF OH (IN cm" 1 ) 


K 

F 

1 

F 

2 

0 

32440. 61 


1 

32474.62 

32474.30 

2 

32542.56 

32541 .98 

3 

32644.22 

32643.45 

4 

32779.^9 

32778.49 

5 

32948.31 

32947’. 05 
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energies for the £ state. In this case F and F are the + and - 

symmetry states of a given N. They are labeled in J according, to. 

J = N + 1/2 for F^ and J = N - 1/2 for F^. The symmetry is labeled 
+ for N = even and - for N = odd in the £ state. The f^ levels are 
+ symmetry for K = even and - for K = odd, while the f ' levels are 
- for K = even and + for K = odd. Figures 1, 2, and 3 show the correct 
labeling of the energy levels. 

D. The Determination of the Transition Probabilities 

A molecule may make a transition from a given energy state to a 
lower one and hence emit a quantum of radiation spontaneously, under 
the influence of another photon, or because of a disturbing collision 
with another particle. If we set aside the third of these for the moment 
the contribution to the intensity of a spectral line from a slab of 
molecules of thickness ds is: 


dl 


nm 

ds 


N he v(A 
n ' 


nm 


+ B I ) 
nm nm 7 


- N he v B I 
m mn nm 


(II-D-1) 


' where the three Einstein coefficients are related by: 

B g B . 

mn n nm 1 

B " g~ ’ A ~ q . _3 ' 

nm m m nm ojrnc v 


( II-D-2) 


As before the frequency is given in units of Inverse wavelength. Clearly 
only the spontaneous coefficient is needed to fully describe the emission. 

This coefficient is .related to the dipole moment matrix element 
of the transition by: 
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(II-D-3) 


n'v'j' 

n'V'j" 


v 3 

3h 


E 

n 'v'j' 
R u __ u T // 

2 


n v J 



g j' 


where the ' denotes the initial state and " the final state, and g 
is the degeneracy of the level. The summation here is over all sub- 

levels which are not distinguished from one another in the spectral line. 
The n refers to the collection of quantum numbers of the state excluding 
the vibration v and the rotation J. 

In the approximation that the wavefunction of the molecular state 
is separable, there is no interaction between the electronic motion, 
the vibration, and the rotation of the molecule, and the matrix element 
is: 


= = (n , ,v / ,J / Q / M / |M e |n /, ,v // ,J // Q ,/ M // ) 

+ <n 'I n"y<v / ,J / n / M / |M n | v",J"q"M ; ') 


(lItD-4) 


where Dirac's notation has been used, and M is that part of the dipole 

moment operator independent of the internuelear distance while M is 

n 

the rest. For n' -JL n" the electronic eigenfunctions are orthogonal and 
the second term above drops out. Thus: 

l 


R n*v # J* = (n%v / ,J / n / M / |M e |n // ,v // ,J / 'Q // M"^ 


(II-D-5) 


Landau has shown that: 


mV' 


R 


n 'v ' J ' 

n'V'j" 



3 H i R J 2 ( n - D - 6 ) 

M 'M" 
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so that if the dipole operator in molecule fixed coordinates is expressed 

as M = M sine + M cos0 then: 
z 71 C 


R z = VM'jsin^JjVM'') 

+ ^n'|Mjn' , ^j'Q / M'|sinecose|j''fi / 'M / ^v / |r 2 |v // ^ . (II-D- 7 ) 


This contains tn only in the form exp(iMo)) so that the second expression 
in each term above will vanish unless M' = M". Similarly, the expression 
of R x ± IRy will lead to the selection rule AM = ±1. If the wave- 
function varies only slowly with the internuclear distance, averaging 
over r allows the Franck-Condon overlap integral to be separated: 




■Sr ■ ( v 'l r2 l v ‘') 

'av 


(ii-d-8) 


and the electric dipole moment: 


R^ V ' = (n'|M|n'\ 
x 'av 


( 11 -D- 9 ) 


The remaining integral is written as: 


1 


T J'fl'M' 
J"Q%" 


= ^j'n'M'jsinecoseljVM") + (j'n'M'Isi^ejjVM^ . (ii-d-io) 


•The remaining selection rules along with the ones deduced above are: 


AM = 0, ± 1 , AS = 0 

AJ = 6, ± 1 , AS = 0 

AQ - 0, ± 1 


(ir-D-n) 
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Finally combining all this gives the matrix element: 



M M 


n 'v ' J r 

2 

n / n // 

2 

v V' 

Vv"j" 


R 

.e 


R . . 

vib 


3 j'j " 


(II-D-12) 


where S t/t// is the line strength and obeys the summation rules: 
J J ' " ‘ 


■ 2J ' + i - E s j'j" - 2J,,+1 • (11-K-13) 


J" 


The transitions allowed by the selection rules between states of the 
OH molecule are labeled 0, P, Q, and R according to the scheme listed in 
Table 6. These labels do not consider the effects of symmetry or hyper- 
fine structure. Table 7 gives the line strengths of the various transi- 
tions of interest calculated from the general formulae of Kovacs, and Table 
8 gives the resulting spontaneous transition probability matrix. In Table 
8 the energy levels are numbered in order of increasing energy refering to 
Figures 1, 2 , and 3: levels I-16 are the II states with V = 0, 17~32 the 
V = 1 levels, 33-48 the V = 2 states and levels are the £ 

-states. 

The above calculations can be easily extended to include hyperfine 
structure if desired by writing: 


[ ( 2J '+1 ) ( 2F // +1)W 2 (J / J // F V' ; II ) ]Sj , 3 „ 


(II-D-14) 


in place of S / „ is the square of the Racah coefficient listed 

<J J 


in Table 5 




TABLE 6 

The 12 transitions allowed are (Dicke and Croswhite 1962 ): 

Name Transition 

0 12 ‘ (K-2, J-l) - (K,J) j 1-2 

P x (K-l, J-l) - (KjJ) ; 1-1 

P 2 (K-l, J-l) - (K,J) ; 2 - 2 

'P 12 • (K-l, J) - (K,J) ; 1-2 

Q x (K,J) - (K,J) ; 1 - 1 

Q 21 (K, J-l) - (K,J) ; 2- 1 

Q P (K,J) - (K,J) ; .2 - 2 

Q 12 (K, J+l) - (K,J) ; 1 - 2 

(K+l, J+l) - (K, J) ; 1-1 
R 21 (K+l, J) - (K,J) ; 2- 1 

' R 2 (K+l, J+l) - (K,J) ; 2 - 2 

S 21 (K+2, J+l) - (K,J) ; 2-1 

• 55 - 

Ref er to Figures 1, 2, and 3* 
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TABLE 7 


The line strengths for the 


^I-^I transitions in OH (dimensionless) 


1 

*1.00 

Ri(K-i) 

Pg(K) 

R 2 (k-i) 

%(K) 

Q 2 (k) 

R 21 (K) 

p 12 ( k + i ) 

V K ) 

« 12 (K+1) 

SgjW 

0 la (K + 2) 

1 

0. 

0. 

2.30 

0.667 

0.031 

o.o 4 o 

0.005 

2 

1.62 

1.29 

i. 4 o 

0.300 

0.045 

0.037 

0.006 . 

3 

2.90 

2.3^ 

0.98 

0.223 

0. 0^9 

0.033 

0.008 

4 

4.04 

3-37 

0.73 

0.188 

0.052 

0.027 

0.007 

5 

5.14 

^•37 

o .57 

0.17 

0.05 

0.02 

0.006 

6 

6.21 

5.38 

0 .46 

0.16 

0.05 

0.02 

0.005 


The line strengths for the JC-TI transitions in OH 


K" 

p i 

P 2 

% 

Q 2 

*1 

R 2 

1 

.587 


.563 

-331 

.162 

.169 

2 

• 795 

.275 

1.063 

.688 

.381 

.356 

3 

1.031 

.538 

1.584 

i.i 4 o 


.556 

4 

1.280 

.807 






;°i2 

P 12 

'^21 

Q 12 

R 21 

S 21 

1 


.331 

. 4 oo' 

.175 

.244 

.050 

2 

.081 

.375 

.319 

.238 

.319 

.075 

3 

.106 

.400 

.331 

.262 

.338 


4 

.113 

• 39 ^ 

.306 





\ 


The line strengths for the 


2 2 

2- 2 transitions in OH 


k" 

R- 

R 21 

R 2 

0 

2.66 

2 .66 

2.66 

1 

1.60 

1.60 

1.60 

2 


1.37 
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TABLE 8 


1 

tnc Nfw sp/imintous i*» assures ^»eix ro* case 300 


CASC NU*3tR 300 54 P4WS 54 CJLU^KS STORAGE mddE 0 PAGE I 


COLUMN 

l 


2 

3 

4 

* 


6 


7 


8 


9 

10 


RF1W 

1 

0.0 


0.0 

0.0 

0.0 

0.0 


0.0 


0.0 


0.0 


'o.o 

O.o 


ROW 

2 

8.545 E -11 

o.c 

o.o 

• 0.0 

0.0 


0.0 


0.0 


0*0 


0.0 

0.0 


0QW 

3 

1.370 F -41 

0.0 

7.0 

0.0 

o.o 


o.o* 


0.0 


0*0 


0.0 

0.0 


ROW 

4 

0.0 


1.3736-01 

1.61 SE -09 

0.0 

0.0 


0.0 


9.0 


0.0 


0.0 

O.A 


ROW 

5 

3.4686-02 

0.0 

0.0 

0.0 

0.0 


0.0 


0.0 


0.0 


0.0 

0*0 


ROW 

6 

0.0 


3.476 P -02 

0.0 

0.0 

1.159 P -0? 

0.0 


0.0 


0.0 


0.0 

0.0 


ROW 

7 

0.0 ■ 


4.410 E -02 

8.9446-03 

0.0 

6.435 F -03 

0.0 


0.0 


0.0 


0.0 

0.0 


ROW 

n 

4.433 F -02 

0.0 

0.0 

8*9446-03 

0.0 


6. 4676-03 

1. 1396-09 

0.0 


0-0 

0.4 


ROW 

7 

0.0 


0.0 

5.179 E -01 

o.o' 

0-0 


0*0 


0.0 


0.0 


0.0 

0.4 


ROW 

10 

0.0 


0.0 

0.0 

5.2116-01 

0.0 


0.0 


0.0 


0.0 


9.5136-09 

0.0 


ROW 

U 

1.737 E -02 

0.0 < 

0.0 

5.579 E -02 

0.0 


0.0 


3-510 E -02 

0*0 


3.1236-03 

0,0 


ROW 

12 

0.0 


1.7416-02 

5.613 E -02 

o.o 

0.0 


0.0 


o.o 


3.5126-02 

0.0 

3*1 05 F -93 

ROW 

13 

0.0 


0.0 

0.0 

0.0 

0.0 


0.0 


0.0 


0*0 


1.255 F OO 

0.0 


ROW 

14 

0.0 


0.0 

0.0 

0.0 

0.0 


0.0 


0.0 


0.0 


0.0 

1.259* 

00 

RC 1 W 

15 

0.0 


0.0 

2.674 E -0 Z 

0.0 

0.0 


0.0 


0.0 


0.0 


0.0 

6.1786-02 

ROW 

16 

0.0 


0.0 

0.0 

2.675 E -02 

0.0 


0.0 


0.0 


0.0 


6*2336-02 

0.0 


ROW 

17 

0.0 


2. 3946 03 

1.5696 03 

0*0 

3.737 E 

01 

0.0 


0.0 


2.745 E 

01 

0*0 

0*0 * 


ROW 

18 

2.39 AF 

03 

0.0 

0.0 

1.569 F 03 

0.0 


3.7366 

01 

2.7466 

01 

0.0 


0-0 

<0*0 


ROW 

19 

1. 202 E 

03 

0.0 

0.0 

9.777 F 02 

0.0 


0*0 


2.3486 

01 

0.0 


1.816 E 03 

0.0 


ROW 

20 

0.0 


1.202 F 03 

9.7806 02 

0.0 

0.0 


0.0 


0.0 


2.3486 

01 

0.0 

1*8166 

03 

ROW 

21 

9.244 E 

01 

0.0 

0.0 

0.0 

0.0 


1.394 F 

03 

2.5516 

03 

0.0 


0.0 

0.0 . 


ROW 

22 

0.0 


9. Z 44 F 01 

3.0 

0.0 

I .345 F 

03 

0.0 


0.0 


7.5516 

03 

0.0 

0.0 


ROW 

23 

0.0 


3.754 E 01 

4. I 45 E 01 

0.0 

1-4116 

03 

0*0 


0.0 


3.117 E 

02 

0.0 

0.0 


ROW 

24 

3*7546 

01 

0.0 

0.0 

4.1456 01 

0.0 


1.4116 

03 

3.1186 

02 

0.0 


0.0 

0.0 


ROW 

25 

0.0 


0.0 , f 

1.6576 03 

0.0 

0.0 


0.0 


0.0 


0*0 


0.0 

S .077 E 

02 

ROW 

26 

0.0 


0.0 

0.0 

1.65 RF 03 

0.0 


0.0 


0.0 


0.0 


5.080 F 02 

0*0 


POW 

27 

4*3446 

00 

0.9 

0.0 

3.674 E 01 

0.0 


-0.0 


1.7556 

03 

0.0 


.2.4536 01 

0.0 


ROW 

28 

0.0 


4.365 F 00 

3.6756 01 

0.0 

0.0 


0.0 


0.0 


1.7556 

03 

0.0 

2.453 F 

01 

ROW 

29 

0.0 


0.0 

0.0 

0.0 

0.0 


0.0 


0.0 


0.0 


1.8376 03 

0.0 


ROW 

30 

0.0 


0.0 

0.0 

0.0 

0.0 


0.0 


0.0 


0.0 


0.0 

1. R 37 F 

03 

ROW 

31 

0.0 


0.0 

4.083 E 00 

0.0 

0.0 


0.0 


0.0 


0.0 


0.0 

3.043 F 

01 

ROW 

32 

0.0 


0.0 

0.0 

4.083? 00 

0.0 


0.0 


0.0 


0.0 


3.0456 01 

0.0 


ROW 

33 

0.0 


1.4706 03 

9.99 &E 02 

0.0 

2.4176 

01 

0.0 


0.0 


1.8246 

01 

0.0 

0*0 


ROW 

34 

1.479 E 

03 

0.0 

0.0 

9.996 F 02 

0.0 


2.4176 

01 

1.8246 

01 

0.0 


0*0 

0.0 


ROW 

35 

7.144 E 

02 

0.0 

0.0 

5.9576 02 

0.0 


0.0 


1.5086 

01 

0*0 


1*1746 03 

0.0 


ROW 

36 

0.0 


7.1456 02 

5.9596 02 

0.0 

0.0 


0.0 


0.0 


1*5086 

01 

0*0 

1-1746 

03 

ROW 

37 

5.3956 

01 

0.0 

0.0 

0.0 

0.0 


8.5216 

02 

1.603 E 

03 

0.0 


0.0 

0*0 


ROW 

30 

0.0 


5.195? 01 

0,0 

0.0 

R .522 F 

02 

0.0 


0.0 


1.6026 

03 

0*0 

0.0 


PCV 

39 

0.0 


?. i 05= Oi 

2.4286 01 

0.0 

8.3166 

02 

0.0 


0.0 


1.6926 

02 

0.0 

0.0 


ROW 

40 

2.1056 

01 

0*9 

0.0 

2.4286 01 

0.0 


8.3166 

02 

1 .8 B 2 C 

oz 

0.0 


0.0 

0.0 


ROW 

41 

0.0 


0.0 

9.695 E 02 

0.0 

0.0 


0.0 


0.0 


0.0 


0.0 

3.1136 

02 

ROW 

4? 

0.0 


0.0 

0.0 

9.6856 02 

o.o 


0.0 


0.0 


0.0 


3.1266 02 

O.o 


ROW 

43 

2.3 106 

00 

0.0 

0.0 

2.085^ 01 

0.0 


0.0 


1.0376 

03 

0.0 


1.4546 01 

0.0 


pn« 

44 

0.0 


2. 399? 00 

7. O 05 E 01 

0.0 

0.0 


0.0 


0.0 


1.0376 

03 

0.0 

1.454 E 

01 

ROW 

45 

0.0 


0.0 

0.0 

0.0 

0.0 


O.o 


0.0 


n.o 


I .09 t 6 03 

0.0 


P^W 

46 

0-0 


0.0 

0.0 

0.0 

0.0 


0.0 


0.0 


0.0 


0.0 

I .092 F 

03 

ROW 

47 

0.9 


0.0 

2*1916 00 

0.0 

0.0 


o-o 


0.0 


0*0 


0*0 

1.7046 

01 

ROW 

48 

0.0 


0.0 

D.O 

2. 1976 00 

0.0 


0.0 


0*0 


0.0 


1.7096 01 

0.4 


ROW 

49 

7.2616 

06 

0.0 

0.0 

0.0 

0.0 


4.3836 

06 

1 .068* 

06 

0.0 


0.0 

0.0 


ROW 

50 

0.0 


5.29 -E 06 

0.0 

0.0 

4.4036 

06 

0.0 


3.636 F 

06 

0.0 


0*0 

0.0 


ROW 

51 

0.0 


7.903* 06 

5 .2906 06 

0.0 

1.1076 

07 

0.0 


0.0 


2.4826 

07 

0.0 

0.0 , 


ROM 

52 

1.6506 

06 

0.0 

0.0 

2 . 140 F 06 

0.0 


1.130 F 

07 

4.567 F 

06 

0.0 


0*0 

0.0 


ROW 

53 

7 .3266 

05 

0.0 

0.0 

5.2326 06 

0.0 


0.0 


1.054 F 

06 

0.0 


4.56 IF 06 

0.0 

i 

ROW 

54 

0.0 


2.2766 05 

1.438 E 06 

0.0 

0.0 


0.0 


0.0 


1.59 36 

06 

0*0 

t .461 F 

06 


3 ^ 



TABLE 8 (cont.) 


CASE 

MJWbFR 300 

54 

ROWS 

54 

COLUMNS 


STORAGE 

«nn* o 


PAGE 2 




COUJ XH 

11 


12 


13 


14 


15 


16 


17 

18 

19 

20 

POW 

I 

0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 

0.0 

0.0 

0.0 

r ru 

2 

0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 

0.0 

0.0 

0.0 

ROW 

3 

0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


.0.0 

0.0 

0.0 

0.0 ' 

ROW 

4 

0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 

0.0 * 

0.0 

0.0 

ROW 

5 

0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 

0.0 

0.0 

0.0 

ROW 

6 

0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 

0.0 

0.0 

o.o 

ROW 

7 

0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 

0.0 

0.0 

0.0 

ROW 

8 

0.0 


0.0 


0.0 


0.0 


0.0 . 


0.0 


0.0 

0.0 

0.0 

0.0 

POW 

9 

0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 

0.0 

0.0 

o.o 

ROW 

10 

0.0 


0.0 


0.0 


0.0 


o.o 


0.0 


o.o 

0.0 

0.0 

0.0 

ROW 

11 

0.0 


0.0 


o.o 


0.0 


0.0 


0.0 


0.0 

o.o 

0.0 

0.0 

POW 

12 

7.2136-10 

0-0 


0.0 


0.0 


0.0 . 


0.0 


0.0 

o.o 

o.o 

0.0 

ROW 

13 

0.0 


0.0 


o.o 


0.0 


0.0 


o.o 


0.0 

o.o 

0.0 

0.0 

**OW 

14 

0.0 • 


0.0 


1.3356-03 

0.0 


D . 0 


o.n 


0.0 

0.0 

0.0 

O.n 

POW 

1$ 

1.0076 

00 

0.0 


1 .1696-03 

0.0 


0.0 


o.o 


0.0 

0.0 

0.0 

0.0 

POW 

16 

0.0 


1.0066 

00 

0.0 


U 1716-03 

2.1626-10 

0.0 


0.0 

0.0 

0.0 

0.0 

ROW 

17 

4.0366 

00 

0.0 


o.o 


0.0 


0-0 


0.0 


0.0 

0.0 

0.0 

0.0 

ROW 

18 

0-0 


4.C35* 

00 

0.0 


0.0 


0.0 


0.0 


5.9026-11 

0.0 

0.0 

0.0 

POW 

19 

0.0 


2.6lO c 

01 

0.0 


0.0 


3.0616 

00 

o.o 


1.1036-01 

0.0 

0.0 

0.0 

ROW 

20 

2 .&ILE 

01 

0.0 


0.0 


0.0 


0.0 


3.0616 

00 

0.0 

l . 108 E -01 

1.05 Z *-09 

0.0 

ROW 

21 

0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


3.1486-02 

0.0 

0.0 

0.0 

ROW ♦ 

22 

0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 

3. 1&8 E -02 

0.0 

0.0 

ROW 

23 

2. 2306 

03 

0.0 


0.0 


0.0 


0.0 


0.0 


0.0 

3.8356-02 

8.2596-03 

0.0 

POW 

26 

0.0 


2.229 F 

03 

0.0 


0.0 


0.0 


0.0 


3.0506-02 

0.0 

0.0 

8.2446-03 

ROW 

25 

1.5386 

01 

0.0 


1 *334*6 

03 

0-0 


0.0 


2.0856 

01 

0.0 

0.0 

4.1526-01 

0.0 

ROW 

26 

0.0 


1-5806 

01 

0.0 


1.8336 

03 

2.088* 

PI 

0.0 


o.c 

0.0 

0.0 

4.1726-01 

ROW 

27 

0.0 


1-539 E 

02 

0.0 


0.0 


P .0646 

03 

0.0 


1.4566-02 

0.0 

0.0 

4.7685-02 

ROW 

20 

1.5406 

02 

0.0 


0.0 


0.0 


0.0 


2.0646 

03 

0.0 

1.4606-02, 

4.8005-02 

0.0 

ROW 

29 

0.0 


0.0 


o.o 


3.0056 

02 

1.0436 

01 

0.0 


0.0 

0.0 

0.0 

0.0 

ROW 

30 

0.0 


0.0 


3-0096 

02 

0.0 


0.0 


1.0436 

01 

O.o 

0.0 

o.n 

0.0 

ROW 

31 

2.084 g 

03 

0.0 


1.4616 

01 

0.0 


0.0 


9.6936 

01 

0.0 

0.0 

2.2186-02 

0.0 

ROW 

32 

0.0 


1.5526 

03 

0.0 


1.4816 

01 

9.6976 

01 

0.0 


0.0 

0.0 

0.0 

2.2176-02 

ROW 

33 

2. 8116 

00 

0.0 


0.0 


0.0 


0-0 


0.0 


0.0 

4.8086 03 

3.156* 03 

0.0 

ROW 

34 

0.0 


2. 8116 

00 

0.0 


0.0 


0.0 


0.0 


4.8096 03 

0.0 

0.0 

3. 15 A * 03 

ROW 

35 

0.0 


1.7516 

01 

9.0 


0.0 


2.1916 

00 

0.0 


2.4206 03 

0.0 

0.0 

1.9496 03 

ROW 

36 

1.7526 

01 

0.0 


0.0 


0.0 


0.0 


2.1916 

00 

0.0 

2.4206 03 

1.950 F 03 

0.0 

ROW 

37 

0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


1.8696 02 

0.0 

0.0 

0.0 

POW 

38 

0.0 


0.0 


0.0 


0.0 


0.0 


.0.0 


0.0 

1.869? 02 

0.0 

o.o 

ROW 

39 

1.3456 

03 

0.0 


’0.0 


0.0 


0.0 


0.0 


0.0 

7.366* 01 

8.5486 01 

0.0 

ROW 

40 

0.0 


1.3446. 

03 

0.0 


0.0 


0.0 


0.0 


7.3676 01 

0.0 ’ 

0.0 

0.5465 01. 

ROW 

41 

1.0096 

01 

0.0 


1.2006 

03 

0.0 


0.0 


1-4096 

01 

0.0 

0.0 

3. 3366 03 

0.0 

ROW 

42 

0.0 


1.009* 

01 

0.0 


1.2006 

03 

1.4096 

01 

0.0 


0.0 

0.0 

0.0 

3.3376 03 

*OW 

43 

0.0 


9.4*7* 

00 

0.0 


0.0 


1.34 BF 

03 

0.0 


8.8266 00 

0.0 

0.0 

7.420* 01 

ROW 

44 

9.4696 

00 

0.0 


0.0 


0.0 


0.0 


1.3486 

03 

0-0 

6.6296 OQ 

7.4 P 36 01 

'0.0 

ROW 

45 

0.0 


0.0 


0.0 


1. 6446 

02 

6.6186 

00 

0.0 


0.0 

0.0 

0.0 

0.0 

ROW 

46 

0.0 


0.0 


1.6456 

02 

0.0 


ft.O 


6.6196 

00 

0.0 

0.0 

,0.0 

o.o 

ROW 

47 

1*1306 

03 

0.0 


8.6076 

oo 

0-0 


0.0 


5.9366 

Oi 

0.0 

0.0 

8. P 37 E 00 

0.0 

ROW 

48 

0.0 


1. 133* 

03 

0.0 


8.8286 

00 

5.953 E 

01 

0.0 


0.0 

0.0 

0.0 

8. 2775 00 

ROW 

49 

0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


5. 3696 05 

0.0 

0.0 

0.0 

ROW 

50 

0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 

3.6206 03 

o.o 

0.0 

ROW 

51 

6.9*46 

06 

0.0 ' 


0.0 


0.0 


0.0 


0.0 


0.0 

5.6576 03 

3.613 F 05 

0.0 

ROW 

52 

0.0 


3.5416 

06 

0.0 


0.0 


0.0 


0.0 


U 129 E 05 

0.0 

0.0 

1.451* 05 

ROW 

53 

0.0 


1.7236 

06 

0.0 


0.0 


4.8756 

05 

0.0 


5.0106 05 

0.0 

0.0 

3.2465 05 

ROW 

54 

5.5336 

05 

0.0 


0.0 


a. o 


0.0 


3.528 E 

06 

0.0 - 

1.5556 03 

9.8356 04 

0*0 


35 



TABLE 8 ( cont . ) 


CASE 

3*19 

54 

ROWS 

54 

CCLUWPIS 


STORAGE 

nnt>c 0 

RAGE 3 





COLUMN 

21 

72 


23 


24 


25 * 

26 

27 

26 

. 

29 

30 

tow 

t 

0*0 

0.0 


0.0 


0.0 


0.0 

0.0 

0.0 

0.0 


0.0 

0.0 

ROW 

2 

0.4 

0.0 


0.0 


0.0 


0.0 

o.o - 

0.0 

0.0 


0.0 

0.0 

ROW 

9 

0.0 

0.4 


0.0 


0-0 


0.0 

0.0 

0.0 

0.0 


0.0 

0.0 

ROW 

4 

0.0 

0.0 


0.0 


0.0 


0.0 

0.0 

0.0 

0.0 


0*0 

0.0 

•OW 

5 

0.0 

0.0 


0.0 


0.0 


0.0 

0.9 

0.0 

0.0 


0.0 

o.n 

*OW 

6 

0.0 

o.o 


0.0 


0.0 


0.0 

0.0 

0.0 

0.0 


0.0 

0.0 

ROW 

7 

0.0 

0.0 


0.0 


0.0 


0.0 

o.o 

0.0 

0.0 


0.0 

0.0 

ROW 

a 

0.0 

0.0 


0.0 


0.0 


0.0 

0.0 

0.0 

0.0 


0.0 

0.0 

ROW 

9 

0.0 

0.0 


3.0 


0.0 


0.0 

0.0 

0.0 

0.0 


0.0 

0.0 

*LHW 

10 

0.0 

o.o 


0.0 


0.0 


0.0 

0.0 

0.0 

0.0 


0.0 

0-0 

ROW 

Il- 

0.0 

0.0 


0.0 


0.0 


0.0 

0.0 

0.0 

0.0 


0.0 

0.0 

ROW 

ls 

0.0 

0.0 


0.0 


0.0 


0.0 

0.0 

0.0 

0.0 


0.0 

0.0 

»ow 

13 

0.0 

0.0 


0.0 


0.0 


0.0 

0.0 

0.0 

0.0 


0.0 

0.0 

ROW 

1* 

0.0 

0.0 


o.o 


0.0 


0.0 

0.0 

0.0 

0.0 


0.0 

0.0 

ROW 

15 

0,0 

0.0 


0.0 


0.0 


0.0 

0.0 

0.0 

0.0 


0.0 

0.0 

ROW 

16 

0.0 

0.0 


0.0 ■ 


0.0 


• 0.0 

0.0 

0.0 

o.o 


0.0 

0.0 

ROW 

17 

0.0 

0.0 


0.0 


0.0 


0.0 

0.0 

0.0 

0.0 


0.0 

0.0 

ROW 

i« 

0.0 

0.0 


0.0 


o.o 


0.0 

0.0 

0.0 

0.0 


0.0 

0.0 

ROW 

19 

0.0 

0.0 


0.0 


O.n 


0.0 

0.9 

0.0 

0.0 


0.0 

o.o 

KOW 

20 

0.0 

0.0 


0.0 


0.0 


0.0 

0.0 

0.0 

0.0 


0.0 

o.n 

ROW 

21 

0.0 

0.0 


0.0 


0.0 


0.0 

o.o 

0.0 

0.0 


0.0 

0.0 

ROW 

22 

5.382E-09 

0.0 


0.0 


0.0 


0.0 

0.0 

0.0 

0.0 


0*0 

0.0 

row 

23 

5. 106E-03 

0.0 


0.0 


0.0 


0.0 

0.0 

0.0 

0.0 


0.0 

0.0 

ROW 

24 

0.0 

5-077E- 

-03 

4.636E- 

10 

0.3 


0.0 

0.0 

0.0 

0.0 


0.0 

o.n 

»OW 

25 

0.0 

0.0 


0.0 


0-0 


0-0 

0.0 * 

0.0 

0.0 


0.0 

0.0 

ROW 

26 

0.0 

0.0 


0.0 


0.0 


4.793F-09 

0.0 

0.0 

0.0 


0.0 

0*0 

stow 

27 

0.0 

0.0 


2.766E- 

*02 

0.0 


Z.905E-O3 

0.0 

0.0 

0.0 


0.0 

0.0 

ROW 

26 

o.o 

o.o 


0.0 . 


2.773E-02 

0.0 

2«6B6E~03 

6.304E-10 

0.0 


0.0 

0.0 

ROW 

29 

0.0 

0.0 


0.0 


0.0 


1.094E 03 

0.0 

0.0 

0.0 


0.0 

0.0 

ROW 

30 

0.0 

o.o 


0.0 


0.0 


0.0 

1 • 01 IE 00 

0.0 

0.0 


2.023F-08 

0.0 

SOW 

31 

0.0 

0.0 


9.0 


0.0 


0.0 

S.214E-02 

7.9«4F-0l 

0.0 


1.171F-03 

0.0 

ROW 

32 

0.0 

0.4 


9.0 


0.0 


5.251E-02 

0.0 

0.0 

7.970F-01 

0,0 

1.096F-03 

ROW 

33 

7.476E 01 

0.0 


0.0 


5.490E 

01. 

0.0 

0.0 

fi.0f.9F 00 

0.0 


0.0 

0.0 

ROW 

34 

0.0 

7. <74= 

01 

5.49ZE 

01 

0.0 


0.0 

0.0 

0.0 

8.C56F 

01 

0.0 

0.0 

ROW 

35 

9.0 

0.0 


4.6V3E 

01 

0.0 


3i646? 03 

0.0 

0.0 

5.211E 

or 

0.0 

0.0 

ROW 

36 

0.0 

0.0 


0.0 


4.644? 

01 

0.0 

3.646F 03 

5.216E 01 

0.0 


n.o 

O.Q 

ROW 

37 

0.0 

2. 791 £ 

03 

4 , 904 E 

03 

0.0 


0.0 

0.0 

0.0 

n.o 


0.0 

0.0 

ROW 

38 

2.741E 03 

0.0 


0.0 


6. 122F 

03 

0.0 

0.0 

0.0 

0.0 


0.0 

0.9 

ROW 

39 

2. T46£ 03 

o.n 


3.0 


6 . 05 l E 

02 

0.0 

0.0 

4.379F 03 

0.0 


0.0 

0.0 

BOW 

40 

0.0 

2.745* 

03 

6.CS4E 

02 

0.0 


0.0 

0.0 

0,0 

4.328E 

03 

o.n 

o.n 

bow" 

41 

0.0 

O.O' 


0.0 


0.0 


0.0 

1.022E 03 

3.1T7F 01 

0.0 


3.6R2E 03 

0.0 - 

ROW 

42 

0.0 

0.0 


0.0 


0.0 


1.022E 03 

0.0 

0.0 

3.176E 

01 

0.0 

3.6A0E 03 

BOW 

43 

0.0 

0.0 


3*5325 

03 

0.0 


4.94 3 F 01 

0.0 

0.0 

3.096F 

02 

0.0 

0.0 

RQW 

4 S 

0.0 

0.0 


3.0 


3.532? 

03 

0.0 

4.942E 01 

3.047E 02 

0.0 


0.0 

0.0 

ROW 

44 

0.0 

0.0 


0.0 


0.0 


3.B33f 03 

0.0 

0.0 

0.0 


0.0 

6.043E 02 

ROW 

46 

0.0 

0.0 


0.0 


0.0 


0.0 

3.834? 03 

0.0 

0.0 


6.048F 02 

0.0 

ROW 

*7 

0.0 

0.0 


0.0 


0.0 


0.0 

6.UIE 01 

3.914F 03 

0.0 


7.977E 01 

0.0 

ROW 

48 

0.0 

0.0 


0.0 


0.0 


6.155F 01 

o.n 

0.0 

3.913E 

03 

0.0 

2.9-JlE ot 

Rnw 

44 

4.0 

2.9895 

O r 

7.2855 

04 

0.0 


0.0 

0.0 

0.0 

0.0 


0.0 

o.n 

90W 

50 

3.008? 05 

0.0 


2.478E 

05 

0.0 


0.0 

0.0 

0.0 

0.0 


0.0 

o.n 

Rnw 

5 1 

7.582E 04 

o.c 


0.0 


1.669? 

OS 

0.0 

0.0 

4 . 74 IE 04 

0.0 


0.0 

0.0 

ROW 

52 

0.0 

7.719* 

04 

3.21 1 E 

05 

0.0 


0.0 

0.0 

0,0 

2.415E 

05 

0.0 

0.0 

ROW 

53 

0.0 

0.0 


7.181E 

04 

0.0 


3.11 IE 05 

0.0 

O.n 

1.176E 

05 

0.0 

0.0 

ROW 

54 

0.0 

0.0 


1.093E 

05 

0.0 


0.0 

1.014E 05 

3.7T5E 05 

0.0 


0.0 

0.0 


36 



TABLE 8 (cont.) 


CAS6 

WJK3FR 300 

54 

•nwS 54 

COLUMNS 

$T0»4C6 

HODS 0 

woe 4 




COLUMN 

31 


32 


33 

34 

35. 

36 

37 

36 

39 

40 

ROW 

I 

0.0 


0.0 


0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

RPW 

2 

0,0 


0.0 


0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0*0 

0.0 

row 

3 

0,0 


0.0 


0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

o.'t 

con 

4 

0-0 


0.0 


0.0 

0.0 

0*0 

o.o 

0.0 

0.0 

0.0 

0.0 

tnw 

5 

0,0' 


0.0 


3.0 

0.0 

0.0 

0.0 

0.0 

0.0 

o.o - 

0.0 

tpW 

6 

0,0 


0.0 


0.0 

0.0 

o;o 

0.0 

0.0 

0.0 

0.0 

0.0 

ROW 

1 

0,0 


0.0 


0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

ROW 

R 

0.0 


0.0 


0.0 

0.0 

0*0 

0.0 

0.0 

0.0 

0.0 

0.0 

1 IIM 

9 

0,0 


0.0 


0.0 

0.0 

0*0 

0.0 

0.0 

0.0 

0.0 

0.0 

ROW 

10 

0.0 


0.0 


0.0 

0.0 

0.0 

0.0 

o.o 

o.o 

0.0 

o.o 

ROW 

11 

0-0' 


0.0 


0.0 

0.0 

o.o 

0.0 

o.o 

0.0 

0.0 

O.o 

ROW 

12 

0,0 


0.0 


0.0 

0.0 

0.0 

0*0 

0.0 

0.0 

0.0 

0.0 

pnv 

13 

0,0 


0.0 


0.0 

0.0 

0.0 

o.o 

0.0 

0.0 

0.0 

o.o 

tow 

1 A 

0-0 


0.0 


0.0 

0.0 

0.0 

0.0 

o.o 

0.0 

0.0 

O.n 

tow 

1* 

0.0 


0.0 


0*0 . 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

tnw 

16 

o.o 


0.0 


0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

tow 

17 

0.0 


0.0 


0.0 

0.0 

0*0 

0.0 

0.0 

0.0 

0.0 

0*0. 

ROW 

10 

0.0 


0.0 


0.0 ' 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

tow 

19 

0.0 


0.0 


0.0 

0.0 

0*0 

0.0 

0.0 

0.0 

0.0 

0.0 

tow 

20 

0.0 


0.0 


0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

tow 

21 

0.0 


0.0 


0*0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

tnw 

22 

0.0 


0.0 


0.0 

o.o 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

tow * 

23 

0.0 


0.0 


0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

ROW 

24 

0,0 


0.0 


0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

O.o 

tow 

25 

0.0 


0.0 - 


0.0 

0.0 

0.0 

0.0 

0.0 

o.o 

0.0 

0.0 

ROW 

26 

0.0 


0.0 


0.0 

o.o 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

tow 

27 

0.0 


0.0 


0.0 

0.0 

0*0 

0.0 

0.0 

0.0 

o.o 

0.0 

tow 

20 

0.0 


0.0 


0.0 

o.o 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

ROW 

29 

0.0 


0.0 


0.0 

0.0 

0*0 

0.0 

0.0 

0.0 

0.0 

o.o 

ROW 

30 

0.0 


0.0 


0.0 

o.o 

0.0 

9*0 

0.0 

0.0 

o.o 

0.0 

tow 

31 

0.0 


0.0 


0.0 

0.0 

0*0 

0.0 

0.0 

0.0 

0.0 

0.0 

ROW 

32 

6. 123E-11 

0.0 


0.0 

o.o 

0*0 

0*0 

0.0 

0.0 

0.0 

0.0 

ROW 

33 

0.0 


0.0 


o.o 

0.0 

0.0 

0.0 

’0.0 

0.0 

0.0 

0*0 

ROW 

34 

0.0 


O.Q 


1.5Z2E-09 

0.0 

0.0 

0*0 

0*0 

0.0 

0.0 

0.0 

ROW 

35 

6.1 10E 

00 

0.0 


8.2326-02 

0*0 

O.o 

0*0 

0.0 

0*0 

0.0 

0.0 

ROW 

36 

0.0 


6. 1 10? 

00 

0.0 

8.2796-02 

4.425F-09 

0.0 

0*0 

O.Q 

0.0 

0.0 

ROW 

37 

0.0 


0.0 


2.5976-OZ 

0*0 

0.0 

0*0 

0.0 

0.0 

0.0 

0*0 

ROW 

38 

0.0 


0.0 


0.0 

2.592F-0? 

0.r\ 

0.0 

1*0266-10 

0.0 

0.0 

0.0 

ROW 

39 

0.0 


0.0 


0*0 

1.4926-02 

1.7476-03 

0.0 

1.1936-04 

0.0 

0.0 

0.0 

ROW 

40 

0.0 


0*0 


1. 5026-02 

0.0 

0.0 

1.7346-03 

0.0 

i.2066-04 

1.346F-10 

0.0 

ROW 

41 

0.0 


4.162E 

^01 

0.0 

0.0 

3.095F-01 

0.0 

0.0 

0.0 

0.0 

0.0 

ROW 

*2 

4.t6*e 

01 

0.0 


0*0 

0.0 

0.0 

3.1036-01 

0.0 

o.o 

0.0 

0.0 

ROW 

43 

4. U06 

03 

0,0 


1*1436-02 

0*0 

0*0 

3*8366-02 

0.0 

0.0 

5.R936-02 

0.0 

ROW 

44 

0*0 


4. 11 IE 

03 

0.0 

1,1436-02 

3.835F-02 

0.0 

0.0 

0.0 

0.0 

5. 0996-02 

ROW 

45 

2.0906 

01 

0.0 


0-0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

o.n 

ROW 

46 

0.0 


2.9916 

Ot 

0.0 

0.0 

0.0 

0.0 

0.0 

0*0 

0.0 

0,0 

tow 

47 

0.0 


1.9416 

02 

0.0 

0.0 

1. 6236-02 

0.0 

0.0 

0.0 

0,0 

o.n 

ROW 

48 

1.9S2F 

02 

0.0 


0.0 

0.0 

0.0 

1.7136-02 

0.0 

0.0 

0.0 

0*0 

ROW 

49 

0.0 


0.0 


1.6826 04 

0.0 

0.0 

0.0 

0,0 

9.314E 04 

0.0 

.0.0 

PCW 

*0 

0.0 


0.0 


0.0 

1.1306 04 

olo 

0.0 

9.379C 03 

0.0 

0.0 

0.0 

ROW 

51 

0.0 


0.0 


0.0 

1.7036 04 

1.1276 OS 

0.0 

2.3566 03 

0.0 

0.0 

0.0 

tow 

52 

0.0 


0.0 


3.5266 03 

0.0 

0.0 

4.5666 03 

0.0 

2. 4066 03 

0.0 

0.0 

ROW 

53 

3 . 32 l £ 

04 

0.0 


1*3716 03 

0.0 

0.0 

1.01 *E 04 

0.0 

0.0 

2.2S1E 03 

0.0 

ROW 

54 

0.0 


2.4045 

03 

0.0 

4.6736 03 

3.0776 03 

0*0 

0.0 

0.0 « 

0.0 

3.4146 03 


37 



TABLE 8 (cont.) 


CASF AUwaER 300 5* C 9WS 54 COl'J"«S. STORAGE HOOF 0 


PAGE 5 


CCLl/rifJ 


«OW 

l 

sow 

'2 

RO W 

3 

ROW 

4 

row 

a 

ROM 

6 

ROW 

7 

p.nw 

a 

ROW 

9 

ROW 

10 

row 

(1 

PQW 

12 

ROW 

13 

Row 

14 

ROW 

IS 

ROW 

16 

ROW 

17 

POM 

IB 

PCW 

19 

ROW 

20 

PCW 

21 

ROW 

22 

ROW 

23 

ROW 

24 

pnw 

25 

ROW 

25 

ROW 

27 

ROW 

23 

ROW 

29 

ROW 

39 

ROW 

31 

Rf)w 

32 

ROW 

3J 

ROW 

34 

ROW 

35 

PCW 

36 

ROW 

37 

POW 

3B 

PCW 

39 

pnw 

40 

ROW 

41 

ROW 

42 

Row 

43 

ROW 

44 

pnw 

45 

POw 

46 

ROW 

47 

ROW 

46 

ROW 

49 

pow 

4C 

POw 

51 

row 

52 

ROW 

S3 

ROW 

54 


41 

42 

43 

44 

45 

0*0 

o.c 

0.0 

0.0 

0.0 

0*0 

o.o 

3.0 

0.0 

0.0 

0.0 

0.0 

3.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.9 

0.9 

0 .0 

0.0 

o.c 

0.0 

0.0 

0.0 

0.0 

0.0- 

0.0 

0.0 

0.0 

0.0 

0.0 

0.9 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

o.o 

0.0 

0.0 

0.0 

0.0 

9.0 

0.0 

0.0 

0.0 

9.0 

o.o 

0.0 

0.0 

o.n 

0.0 

0.0 

0.0 

0.0 

0.0 

o.o 

0.0 

0.0 

o.o 

o.o 

0.0 

o.o’ 

O.n 

0.0 

0.0 

0.0 

o.o 

o.o 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

o.o 

0.0 

o.o 

o.o 

0.0 

0.0 

0.0 

0.0 

0.0 

o.o 

0.0 

0.0 

0.0 

0.0 

0.0 

o.o 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

9.0 

0.0 

0.0 

o.o 

0.0 

0.0 

0.0 

0.0 

0.0 

o.o 

3. a 

0.0 

o.c 

o.o 

0.0 

*0.0 

0.0 

0.0 

0.0 

0.0 

o.o 

0.0 

0.0 

0.0 

o.o 

3.0 

0.0 

0.0 

o.o 

0.0 

3.0 ‘ 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

o.o 

0.0 

0*9 

0.0 

o.o 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

c.o 

o.o 

o.o 

0.9 

3.0 

0.9 

0.0 

0.0 

o.c 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

o.o 

c.o 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.9 * 

0.0 

0.0 

'0,0 

0.0 . 

0.0 

0.0 

0.0 

0.0 

0.0 

o.o 

0.0 

o.o 

0.0 

4.997E-C9 

o.c 

3.0 

0.9 

0.9 

2.537E-C3 

0.0 

0.0 

0.0 

o.o 

0.0 

2.S18F-03 

1.394E-1D 

0.0 

0.0 

7.453E-01 

o.o 

o.o 

0.0 

0.0 

0.0 

7.507P-0! 

0.0 

0.0 

1.B96F-05 

0.0 

3.7728-02 

5.037E-91 

0.0 

7.737E-04 

4.133F-02 

0.0 

3.0 

5.B66C-01 

0.0 

'2.767E 03 

o.c 

0.0 

0.0 

0.0 

0.0 

7.726F 03 

0.0 

O.o 

0.0 

0.0 

S.259F 03 

0.0 

0.0 

1.464F 03 

9.697E 03 

0.0 

0.0 

0.0 

0.0 

9.69BE 03 

0.0 

0.0 

0.0 

0.0 

0.0 

3.164P 03 

I.176E 04 

0.0 

0.0 


46 

47 

48 

49 

30 

0.0 

0.0 

0.9 

0.0 

9.0 

0.0 

0.0 

0.0 

0.0 

0.0 

o.o 

0.0 

0.0 

0.0 

0.0 

0.0 

o.o 

0.0 

0.0 

O.Q 

0.0 

0.0 

0.0 • 

0.0 

9.9 

0.0 

0.0 

0.0 

0.0 

0.0 

O.'O 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.9 

0.0 

0.0 

0.0 

0.0 

0.0 

o.o 

o.n 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.9 

0.0 

o.n 

0.0 

0.0 

0.0 

0.0 

0-0 

0.0 

0.0 

0.0 

0.0 

o.n 

9.0 

0.0 

0.0 

0.0 

0.0 

0.0 

o.o 

o.o 

o.o 

0.0 

*0.0 

0.0 

o.o 

0.0 

o.c 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.9 

0.0 

o.o 

0.0 

0.0 

0.0 

0-0 

0.0 

0.0 

o.n 

0.0 

* 0-0 

0.0 

o.o 

0-0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

o.o 

0.0 

0.0 

0.0 

0.0 

o.n 

O.n 

0.0 

0.0 

0.0 

0.0 

0.0 

0.9 

0.0 

0.0 

n.o 

0.0 

9.0 

0.0 

o.o 

0.0 

o.o 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

o.o 

0.0 

0.0 

0.0 

0.0 

0.0 

o.o 

0.0 

o.o 

0.0 

0.0 

o.o 

0.0 

0.0 

0.0 

0.0 

o.o 

0.9 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

o.n 

0.0 

0.0 

0.0 

0.9 

0.0 

0.0 

O.Q 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

o.n 

0.0 

o.o 

0,0 

o.o 

0.0 

0.0 

0.0 

0.0 

0.0 

0.9 

0.0 

0.0 

o.o 

0.0 

c.o 

0.0 

0.0 

0.0 

o.o 

0.0 

o.o 

0.0 

0.0 

0.0 

0.0 

3.0 

O.o 

0.9 

O.o 

0.9 

0.0 

o.o 

0.0 

0.0 ' 

0,0 

0.0 

0,0 

o.o 

o.o 

0,9 

0.0 

0.0 

0,3 

0.0 

0.0 

9.653E-04 

3.5446-06 

0.0 

0,0 

o.n 

0.0 

0.0 

0.9 

0.0 

n.o 

0.0 

0.0 

0.0 

2.196E-02 

o.n 

0.0 

0.0 

0.0 

t. 018-02 

0.0 

7.617F 03 

0.0 

0.0 

o.o 

8.910E-02 

3.6478 03 

1.034F* 03 

0.0 

o.n 

o.n 

0.0 

0.0 

7.4&6E 03 

0.0 

0.0 
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TABLE 8 ( cont . ) 


CASF fcUMRgR 300 5^ *P«S 5*. COLO“WS 

ccu^ N 5l 52 * 3 54 


ROW 

1 

0.0 

pnw 

2 

0.0 

ROW 

3 

0.0 

‘ ROW 

4 

0.0 

ROW 

3 

0.0 

ROW 

6 

0.0 

ROW 

7 

0.0 

ROW 

a 

0.0 

ROW 

9 

0.0 

ROW 

to 

.0.0 

ROW 

u 

0.0 

ROW 

12 

0.0 

RCW 

13 

0.0 

ROW 

14 

0.0 

J»OW 

15 

0.0 

pnw 

16 

0.0 

'row 

17 

0.0 

ROW 

te 

0.0 

ROW 

19 

0.0 

ROW 

20 

0.0 

ROW 

21 

o.o 

ROW 

22 

0.0 

ROW 

21 

0.0 

ROW 

24 

0.0 

ROW 

2* 

0.0 

ROW 

26 

0.0 

Pnw 

27 

0.0 

ROW 

28 

0.0 

ROW 

29 

0.0 

ROW 

30 

0.0 

PQW 

31 

0.0 

ROW 

32 

0.0 * 

ROW 

33 

0.0 

ROW 

34 

0*0 

ROW 

35 

0.0 

ROW 

36 

0.0 

ROW 

37 

0.0 

POW 

38 

0.0 

ROW 

39 

0.0 

ROW 

40 

0.0 

ROW 

41 

0.0 

POW 

42 

0.0 

ROW 

43 

0.0 

ROW 

4 * 

0.0 

»0W 

45 

0.0 

POW 

46 

0.0 

ROW 

47 

0.0 

ROW 

48 

0.0 

ROW 

49 

0.0 

ROW 

50 

0.0 

ROW 

51 

0.0 

ROW 

52 

1.601E-01 

ROW 

53 

i.09*f-0l 

POW 

54 

0.0 


0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

o.o 

‘0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

o.o 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0,0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

o.o 

0.0 

0.5 

o.o 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

o.o 

0.0 

. 0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0-0 

0-0 

0.0 

0.0 

0.0 

o.o 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

o.o 

0.0 , 

' 0.0 

0.0 

0.0 

3.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.3 

n.o 

o.o 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

o.c 

0.0 

0.0 

3.0 

0.0 

0*0 

o.o 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

o.o 

0.0 

0.0 

0.0 

0.0 

0.0 

o.o 

0.0 

0*0 

0.0 

0.0 

0.0 

o.o 

0.0 

3.0 

0.0 

0.0 

0.0 

0.0 

0*0 

3.0 

o.o 

0.0 

o.o * 

0.0 

0.0 

3.0 

o.o 

0.0 

0.0 

o.o 

c.o 

3.0 

0.0 

0*3 

0.0 

0.0 

0.0 

0.0 

o.o 

0.0 

0.0 

0.0 

3.612**-0 l 

S .077 E -0: 

0.0 
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E . The Population and Transfer Equations 


The energy state of an OH molecule can change by one of four possible 

mechanisms. First the molecule may spontaneously emit a photon and change 

to a state lower by the energy of the photon. Second it may encounter a 

photon and be induced to emit or absorb a photon of the same energy and 

change to a different state. Third xt may encounter a particle of finite 

mass and be induced to change its energy state without emitting or 

absorbing a photon. Finally, it may change states under the influence 

of a particle and emit a photon. The rate of change of the population 
th 

of the i level then is.: 


where 


dN. 

= 2 W..N. - 
dt . Ji J 

3 


N. 

1 


SW i3 

J 


(II-E-1) 


w ik 3 < A ik + V,, + °ik + Vi “ d w ii = 0 

lk Ik 

and' A n , , B,, , Cl, , and D, , are the rate coefficients for the four 

lk lk 7 lk' lk 

mechanisms above in order, (j) is the line shape factor, which will 

v lk 

be discussed later. For our purposes the A and D factors may be combined 
as A^ = A^ + Degeneracy may also be included by writing: 


* 

B. . = g.B. . = g.B.. 
1J i ij j J-i 


A. . = g. A. 7 . 

3. 1J 


* N i 
N. = — 

1 S i 


(II-E-3) 


This results in an equation similar to (11-E-l) with the N’s and W' s 

'X' 

replaced by the corresponding N ’ s and W 's, and the left hand side 

multiplied by the degeneracy factor g_^ . 

Since only two specific levels of a complicated energy structure 

* 

are of interest in a line transition, the rate equations can be written 
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specifically separating those two levels; the subscript u for the upper 
and for the lower of them the two rate equations are: 


dt 

dN £ 

dt 


■Wvi - ' Vie + V-Vi - Wt ’ 

(II-E-4) 

-V-eu^i + VWA + Vud’e + ¥ n - b 4uH s - , 


where 




R u " ^ W ju\N-N, 

'tu 



r = s w . - w o , 

A u ■ . uj u-L ’ 

J 

T i - s % “ v * 

J 

N £u=V N u * 


(II-E-5) 


N = J H. . 

j J 


The subscript of (j) distinguishes the evaluation of <j) v for emission 

ij 

e, or absorption a, or one of the relaxation processes R or p. 

Integrating the rate equations (II-E-4) over frequency, averaging 
over direction, and applying the transformation (lI-E-3) the two rate . 
equations of interest become: ‘ 


4i 



(II-E-6) 


E. 


dN 
u 

u dt 


* *, « 


* * 




* * 


/ A A » A /k A A . A A « A A 

= \ )B J - N u A + R u (N -N^ - , 


dN 


* 


* 


* 


* * 




* * 


C , A A , A A A A , A A i A A 

s e -ST - -(Nt - \ )B J - N u A 4- Bf, (K -K^j) - T t \ , 


•where I, 3, and (j) are related by: 


J OO -CO c 

<t> v dv - 1 > J 4v I v CIV = 1 > J 


I^dv = lAv , 


I 


TSJ IdC2 = J * 

kn 

(ri-E-7) 


The above equations cannot be used to determine the actual 

populations of the levels considered because the relaxation term Ih 

depends upon the populations of the other levels of the molecules 

energy structure. A more workable approach to be used here is to solve 

the simultaneous equations (11-E-l) along with a normality condition such 

as E N. = 1 for the populations -and hence determine the values of the 
3 3 

R's by: 


R = 
u 


= 


N C —NpRp 

U U -V -LU 


N-N 


■in 


NnCn-N R n 
■v -t, u U'C 


N-N 




(II-E-8) 


with 


C, s s 1 y k s= u,t 
k ^ km 3 3 

m 

which can then be used in the transfer equation solution. 

The equation of transfer for the radiation from the two states has 
the form: 
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(II-E-9) 


dlv hy. . 

ST = T? HW, + N uU% - W,U ' 


In the following we will assume that (j) has the same form for emission 
and absorption. Differences between the two functions have not been 
mentioned in the literature, but could arise from the line narrowing 
characteristic of unsaturated maser amplification or other processes.. 

By integrating over frequency and combining with the steady state solu- 
tion of the population equations (II-E-6), the transfer equation may be 
written in dimensionless form, after Goldreich and Keeley (1972): 

—• = ^ l- 1 / .- 2 ) , + s (II-E-10) 

ds ' ' 


where 


3 


BI 

* 

A 



R = 


# 

R + 
u 


* 

% 


R 


r P R +r Rp 

■i, u U -j. 

R +Rp 
u £ 


a = 


* 


RA 




-r 


Rp 

u 



Qg 

l-a ’ 



4jrAv l 

* AN * 
B hv 0 o 


* * 
N +Nn 

,. T u -t 

AN o " 2 S 


(iI-E-ll) 


In the above relations, 3 and 5 are just the dimensionless intensity 

and average intensity, R the total rate of input to the two levels, T a 

the average relaxation rate, and a and p are the pumping and decay 

rates of the pair of levels. L is the unsaturated growth length and 

AN is the unsaturated population difference. S is a term including 
o 

the spontaneous emission, and is related to the excitation temperature 
of the levels in the following way. Anticipating a population inversion, 
the excitation temperature may be written as: 
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r u / hv ut\ 

r exp \sd 


( II-E-12) 


h AH o 

For ^ « 1 and — — « 1 the population difference gives: 


AN . kT 

° = -L and S =* ^rr~ » 1 . 


M 2S 


hv 


ul 


(II-E-13) 


The equation of transfer (lI-E-10) can be solved exactly only in 
one dimensional geometries. Analytical solutions for higher dimensioned 
situations can be found only for unsaturated conditions, where 5 « {3. 
The following chapter will discuss the solution of '(II-E-IO) for 
spherical geometries. 
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CHAPTER III 

THE SOLUTION TO THE TRANSFER EQUATION 

To analyze the population inversion in the states of any molecule, 
the rate equations expressing the total rates of entry to and exit from 
each energy level must be solved. These equations must include both 
collisional and radiative processes and depend upon the properties of 
the emission region. Since the radiation intensities will change through 
the region, so will the populations of the energy levels. A study of 
population ratios of close pairs of levels must consider simultaneously 
the solution of the equations of statistical equilibrium and of radia- 
tive transfer even for the simple case of time independent emission. For 
the case of the OH molecule there are over one hundred energy levels 
relevant to the study of the maser emission; this presents a formidable 
computational problem to solve the emission equations by this brute 
force method. 

A more usable approach is to consider only the needed pairs of 
levels for the spectral lines of interest. Separate analyses for differ- 
ent lines can then be combined in various ways to provide a qualitative 
understanding of the interactions of the multiple rate and transfer 
equations. Quantitative calculations of the stimulated emission are of 
doubtful value with this approach because of the complexities of the 
interactions. Chapter V describes the efforts ,of other authors in 
solving the transfer equations along these lines. Our approach will be 
to solve the complete rate equations coupled with the transfer equation 
for one maser line. 

If stimulated emission is not very important in the other energy 





levels, and if the induced transitions in the pair of levels considered 
does not affect the populations of the other levels, this simplification 
can be completely satisfactory. Our solution will be useful if the 
transition considered is the dominant transfer path out of an energy 
level, and can be applied quantitatively to more complex situations 

by iterating the rate equations with the changing intensities. The effects 

# 

of the interactions are explored in Chapters IV and V. In this chapter 
we will explore the coupled solution of the rate and transfer equations 
for the two levels of an inverted population ratio using equations 
II-E-8, -10, and -11 for a simple model gas cloud. The treatment of 
Goldreich and Keeley (1972) will be extended to generalized degenerate 
energy levels and to all degrees of saturation, correcting some signifi- 
cant omissions in their derivations, and removing some simplifications. 
Their results for partially saturated masers omit several terms and 
use an incorrect approximation of the transfer solution for small core 
sizes, which leads them to the false conclusion of a constant minimum 
core size. They also ignore the effects of the background radiation, 
which we will show to be important for microwave masers in determining 
the saturation state. The same omissions mislead them to find a core 
of constant intensity in fully saturated masers. Further, their 
expressions over simplify the effects of the population inversion ratio 
on the maser characteristics. In the following sections we will derive 
the relations among the maser characteristics- avoiding the errors and 
simplifications, and in Chapter IV show that the solutions presented give 
an understandable and interpretable picture of the nature of the maser 
process. 
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A. Introduction 


One of the simplest models for a cosmic maser is a homogeneous 
spherical cloud of molecules in which inverted population is created at 
a uniform rate. Even in such a simple model it is not necessary that 
the size of the cloud be the observed size of the masing region. Rather, 
for moderate degrees of saturation, tfie observed source could be sub- 
stantially smaller than the actual size of the region. Additionally, 
even if the inversion is created uniformly, it will not be uniform 

in the cloud. If 10° K isotropic external radiation is incident on 

12 

a sphere with a gain of 10 . through a diameter, the radiation density 

7 1*3 

will be 10 at the center and roughly 10 times the fractional solid 
angle of the outward radiation at the edges of the cloud. If the 
angular distribution of the radiation is more than a few milliradians 
wide the radiation density and hence the saturation will be largest , 
at the surface. This phenomenon has been observed in other one dimen- 
sional geometries .in laboratory lasers, and has been suggested by many 
authors to occur in astrophysical masers as well. The apparent reduc- 
tion in source diameter and the increase of saturation toward the outer 
edges of a spherical region are shown to be related in the following 
sections . 

In the following three sections we will solve the transfer equation 
for the limiting cases of almost no saturation, partial saturation with 
a central core of unsaturated amplification, and total saturation. The 
results are presented graphically in Chapter IV. 

The derivations here are restricted in several ways. Since the 
transfer equation explicitly contains the solution of the two-level 
population equations given in the previous chapter, the solutions pre- 
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sented here entrain all the eariier 'restrictions . One significant 
restriction is that the line shapes for emission and absorption be the 
same; if the line of interest is an incomplete overlap of several lines, 
this may not be true. Polarization, scattering, and Zeeman effects will 

l 

be neglected and with them, any cross relaxation effects. Inclusion 
of the hyperfine structure and. the Zeeman structure in the rate equations 
would lift some of these conditions, but the anticipated width for the- 
infrared lines exceed the energy gaps even with line narrowing by maser 
amplification . 

• The equation of transfer from Chapter II is: 


. §-(&)»*■♦!&) • ■ <«> 

The most general solution to this equation, assuming that the 5 term 
is only a function of the linear dimension and not of the intensity, is: 


3 ( s ) = 3 ( s 0 ) 



P 

P+5 





ds 


■] 


dy 


+ 



(III-A-2) 


If 5 is a function of intensity then this is just an integral form 
of (III-A-1). 

» 

B. Solution of the Transfer Equation for No Saturation 

In solving the transfer equation in regions of no saturation the $ 
term from the population equations is set to zero,- so that the stimulated 
emission has no effect on the populations of the maser levels. This will 
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occur only if the pumping rate is high enough. If $ is set to zero in 
the above equation the terms may be integrated .to find: 


3(s) = 5(s o )exp[s-S Q ] + (S + 1/2) [exp[ s-s o ]-l} , (III-B-1) 


where s is the path length in the medium (refer to Eq, II-E-ll) . 

For a spherical geometry a ray path through the region along the 
line of sight will make an angle Q with a radius vector. If p = cos0, 
the ray path length is s q = -Rp to s = +Rp with R the total radius 
of the region. 



Using these limits, the intensity in the positive direction from 
(III-B-1) is a function of the total radius R and the offset from the 
region’s center and is given by: 


3(R,p) = .l(R^-l)exp[R(p+jp|)] + ( S+l/2) £exp[ R( p+ 1 p | ) ] -1 } (lII-B-2) 

where 3(R,-l) is the isotropic background radiation intensity of the 
appropriate frequency. 

Hie average intensity at a point r in the spherical cloud is 
given by II-E- 7 , with dQ = sin0d0d<[) = dpd<|>, or: 


, +1 2it • +1 

= \) ~ I f l(*/n)dpdj) = 5 I I ( r ^) d l J 
^1 ^0 -1 


(III-B-3) 
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The resulting average intensity for totally unsaturated regions is: 



(III-B-4) 


where is the background flux 3(R, -l) incident in the rear of the 

cloud. 


Equations (lII-B-2) and (lII-B-4) will be valid as long as the 
condition <j. « p remains valid'. The maximum size of such a region 
may be established by replacing the inequality with an equality, so that 
the stimulated rate of population transfer is equal to the average 
natural decay: 


BJ = | (A*T a ) 


(III-B-5) 


In the dimensionless forms of the preceding chapter this condition 
becomes: 

$ . (III-B-6) 

Setting the intensity given by (III-B-^) equal to £ gives the radius 
which is the upper limit for the size of a totally unsaturated re- 
gion from the relationship: \ 


1 _ 1 
a ~ 2 





(III-B-7) 


where * l/a is given by the pumping mechanism. 

The dimensionless parameters in the preceding equations may be 
related to the physical conditions of the maser region with the following 
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equations. The physical size is related to the parameters of the region 
by: 


* 

C D = D 


x ^Ao-^vw \i\ 


[l sec ~ 1 j\ 100 ^f \{&/vJ\ 1 cm 3 7 \ N 0 H / \i.5i x 10 12 ' 

where D is the physical diameter in cm, and 


emy 


r = 


T mp-J 

| i+3°k 

for 

• X = 100p 

v k 

1 0 . 08 °K 

for 

X = l8 cm 


1 ( l+r\ 


T 

B 

U 

si 

I* 3 

z v-*l 

} 

1^ 

11 

C") 


* 


[ 


i - 1 

a 


•= i +. 


a j 1-r. 
A /ll+r 


(III-B-8) 


The saturation boundary and the intensity then are dependent on three 

* ' 

parameters: the scaled size of the region C D, the population ratio r, , 

and the background intensity 5 = T /T . The boundary relation becomes: 

o o v 



(III-B-9) 


and the intensity at the center of the observed disk is 



T 

v 



e 


2R 


+ 



( III-B-10) 


The distribution of intensity across the disc given by equation 

(III-B-2) is essentially gaussian. The half width of this is related to 

1/2 

the apparent source size, and is proportional to R . Using the scale 

’X' 

length (C D), the ratio of real to apparent size is: 


D * 

real _ C D 

D app (C*D) 1/2 


(C*D) 1/2 


(1II-B-11) 
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C. Partially Saturated Spherical Clouds 


If the radiation density at the outer edges of the spherical cloud 
of OH exceeds p, or equivalently if the radius of the cloud exceeds 
as given in the preceding section, at least part of the cloud will 
be saturated. Saturation is characterized by a high radiation density 
and hence high stimulated transition rates, which depopulates the upper 
of the two states involved in the maser, and in the extreme forces 
the two populations to become equal. The transfer equation is charac- 

4 

terized by small gain so that: 

xio C( s+1 /2)(e X -l)] = (S+l/2)x (HI-O-1) 

and the gain is proportional to the path length. Logically, higher 
pumping rates in the maser mechanism will reduce the degree of satura- 
tion if all other parameters are constant, and lower rates will increase 
saturation. if the pump rate is a constant, then larger values for the 
cloud radius will give higher radiation densities and greater saturation. 

We can model this saturation behavior in a spherical cloud as a 
core of radius a which is unsaturated surrounded by a shell of saturated 
population. The intensity in the shell is directed mainly outward with 
a mean value of ^(a) = p at the boundary between the core and the 
saturated shell. The gain outward of the boundary is proportional to 
the path length, so that the transfer equation may be solved using: 

*(*) - + p ' (III-C-2) 

from Goldreich and Keeley ( 1972 ), where the first term is the outward 
component of the boundary intensity times the linear gain, the second 


52 




term is the correction for measuring r from the center of the sphere, 
and the third term is the additional contribution of the core in the 
spherical geometry. Lang and Bender (1973) consider a similar approxi- 
mation for the case of complete saturation as a first estimate and solve 
the transfer equation and coupled rate equations iteratively using 
numerical methods. Twenty iterations change the above expressions by 

less than twenty percent and we conclude that (lII-C-2) is a reasonable 

/ 

approximation for the gain behavior in saturated regions. 

The characteristic integral to be evaluated in the transfer equa- 
tion solution (III-A-2) is: 


pds f 3r 2 dr 

^ 3 r 3 +3r 2 +a 2 (3-a) 


(III-C-3) 


The third order polynomial in the denominator of this integral makes 
exact evaluation very difficult, but the problem can be highly simplified 
by a few reasonable approximations. If the gain rate in the saturated 
region is large .enough then (3 « '^ will hold and the integral becomes: 


(3ds f gds f 3r 2 dr 

W-jT-J r 3 + ,2 (3 .,) 


(iii-c-4) 


This will be shown to become unworkable for small values of the core 
radius a ( see page 53) ; for small values of a the last term in the 
denominator of (lII-C-3) is negligible, and the integral becomes: 


Jp? 3 ds ifife dr • 

These two integrals are easily evaluated; the first is: 


(III-C-5) 
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(in-c-6) 



3r 2 dr 

r 3 +a 2 (3~a) 



s 3 +a 2 (3-a) 

s o 3 +a 2 (3 _ a ) 


and the second is: 



(III-C-7) 


Goldreich and Keeley incorrectly apply ( III-C-6) to the limit of small 
a , hence finding a discontinuity in the core size as complete saturation 
is approached. 

Relationships for the core size and the total intensity may be 
found by integrating along a line of sight through the center of the 
cloud using equation (III-A-2). At the rear core boundary: 


3(a,-l) = 5(B,-l) 


R 3 +3a 2 -a 3 

. s 

r 4 4 ,1 

R , „ 2„ 3„ _ 3, 4 

_ 2 

„ 2 

■jj" — 1“ 3* a R' — • -3a 4" a 

L 3a 

3 a 

»»«. _ 


1 3 3 

+ — [R 


(in-c-8) 


for the approximation (lII-C-4). If the s gain in the rear side of the 
cloud is not too large the relation 3(a,~p) = 3(a ? -i) holds, 
hence 3(a^+p) is given by: 


5(a,+fji) = 3(a,-l)e 2a ^ + (S+l/2)(e 2afI -l) 


and 


5(a) = 


3(a,-l) 

f 2a 
e -1 

2 

2a 


+ 1 


§H) 


2a 

e -1 
2a 


(III-C-9) 


-1 


Combining (III-C-8) and (III-C-9) gives the relationship between the 
core size, the total size, and the pumping rate, with 3(a) = g: 
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The integration path is completed by applying (III-C-8) to find the 
surface intensity on the line of sight: 


3(R,1) = fJ(R,-l) 


R J +3a -a' 


L 3 a 


32 

2a j 

^R 3 -a 3 \ 

[ / 

e + 1 

^a 2 ) 

1 + | 
- ' 


R 3 +3a 2 -a 3N 

3a 2 


2a 


+ ~~ |(R 3 +3a 2 -a 3 )(e 2a -l) + 

3 a 


r 4 4 

1 R 2 3 a 3 4 

+ 3a R-a°R - - 3a°+a 


i + 1 y 

3a" 


e 2a -l\ / R 3 -f 3 a^-a 3 ^ 

, 2 /\ 3a 2 


(Ill-C-Il) 


It can be seen that in the limit R = a the equations ( III-C-lo) 
and (lII-C-11) become equations (III-B- 9 ) and (III-B-2) respectively. 
For small values of a (III~C-10) looses meaning as the only solution 
for a = 0 is the trivial R = 0 for all a and- (HI-C-ll) becomes 
useless. Furthermore, equation (Ill-C-Il) is cumbersome to use in 
finding the brightness temperature . 

The second approximation ( III-C- 5 ) gives an indication of the be- 
havior of the cloud as the saturation becomes complete. Again we inte- 
grate along the line of sight on a diameter of the cloud . For the back 
half of the saturated shell the intensity is: 
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_ 4 r 3+R~| 5 . s[/3+r\ 4 ..1^, , ^ ir/3+R\ 3 " 
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The intensity crossing the core is again given by equations (III-C- 9 ) 
so that the intensity at the core boundary is: 


‘ 2 a 
e -I 


(III-C-13) 


Setting $(a) = p again gives the relationship between the core 
size, the total size, and the pumping rate as: 



The peak intensity is given at (R,+l) by: 



(III-C-I5) 

and once more these formulae reduce to the unsaturated case when R = a. 

We can now determine the limit for the cloud to retain an unsaturated 
core by setting R = R^ at a = 0 in equation (HIyC-14) above-. The 


56 



is; 


resulting relationship for R,. 




- i 


(m-c-16) 


Again, equations (III-C-14) and (III-C-I 5 ) become equations (III-B-7) 
and (lII-B-2) when R = a. Furthermore, the solution for the first 
approximation matches the solution for the second approximation at 
R = a, so that for small values of R the second approximation will 
best represent the maser emission. This will also apply for large values 
of R if the saturation is almost complete and equations .(III-C-14) arid 
(II I-C-15) will describe the emission. Otherwise the first approximation 
is the appropriate choice and equations (III-OIO) and (Ill-Oil) should 
be used. The core size is given by (111-015) or (III-OIO) respectively. 

The distribution of intensity is no longer gaussian if the cloud 
is partially saturated. In the region of the disk near the center 
the gain through the saturated shell is nearly constant whether or 
not the line of sight passes through the core. However, the core contri- 
butes an extra (exp(2a)~l) to the gain so that even for small core 
^sizes the core dominates the apparent size of the cloud. Since the 
core's contribution is roughly a gaussian function, the apparent size 
is given by: 


D real _ C* D 
D app ( 2C* L ) 1/2 


(III-C- 17 ) 


where L is the physical size of the unsaturated core. Thus, clouds 

l 

that are appreciably saturated can appear much smaller than their 
physical size. 
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D. Completely Saturated Masers 

If R is greater than the value ' given by equation (lII-C-16) 
the region is over saturated. The intensity at the center of the cloud 
exceeds the boundary value of 5(0) = (3 by a factor K > 1. The average 
intensity $(r) in the transfer equation may be written: 

$( r ) = K . (iII-D-l) 

To find the solution of the transfer equation in this case we apply 
(iII-D-l) to the integral (lII-A-2), and follow the technique of Section 
C. The over-saturation parameter K is related to the other variables 
by setting 3(0) = |3K, giving: 



The peak brightness temperature is given at the center of the disk by: 



With the gain in the maser roughly proportional to the path length, 
the saturated maser has a Lorentzian intensity distribution and the half- 
intensity points are at about 0“-7 R* Then the apparent size is related by: 

p- -— ■ ■ =1.4 . (III-D-4) 

app 
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Lang and Bender (l9T3) give a value of about 3.8 for the above ratio by 
considering in detail the intensity distribution across the disc. 


E. Estimating the Emitted Intensity 

Equations (lII-B~lO), (lII-C- 15 ), and (lII-D-3) give the maximum 
brightness temperature of the maser emission, at the center of the 
observed disc. This could be converted into an intensity if the size 
of the disc were known in sterradians with the standard formula, but 
since the OH cloud will probably be unresolved in the receiving 
antenna, the effective brightness temperature of the entire disc is 
needed. The effective brightness is essentially the average brightness 
temperature over the physical disc, and for unsaturated masers is given 
by: 


< t b> - t b 


B I max 2 c"d( 1 -r) 
or for fully saturated masers by: 


(iII-E-l) 


< t b)' t b 


• max 


CT D( 1 -r) 


(III-E-2) 


Partially saturated masers with a core of size L have an average 
brightness temperature of approximately: ^ 


<v> - t b 4 

max L 


2C*L( 1-r) 


* , 
C D( 


H 

(l-r)J 


(III-E-3) 


Using these the intensity of the maser emission can be estimated from 
equations (iV-D-l) and (lV-D-2). 
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CHAPTER IV 

A GRAPHICAL SOLUTION OF THE TRANSFER EQUATION 
A, Introduction . 

The theoretical analysis of molecular line emission from inter- 
stellar sources involves the simultaneous solution of the equations for 
line transfer and population equilibrium. Most -of the emission of 
astrophysical interest is very simple and any number of simplifying 
assumptions can be introduced in the solution of the coupled rate and 
transfer equations without loosing physical significance. The key to 
these simplifications is the assumption of local thermodynamic equili- 
brium; remarkable agreement between theory and observation has been' 
achieved under this assumption for a wide variety of astrophysical 
sources. For a few radio emission sources involving fine structure 
emission from OH and other molecules, this assumption falls far short of 
adequate. Brightness temperatures exceeding 10^ °K make any LTE 
considerations untennable for these sources; the only currently acceptable 
explanation for the intense radiation is amplification by stimulated 
emission. 

This maser action makes theoretical analysis very difficult because 

v 

the population and transfer equations become interdependent and very , 
nonlinear. Various attempts at the solutions have involved ignoring 
the rate equations while solving the transfer equations for stimulated 
emission, approximating a transfer equation in population solutions, and 
solving a coupled two level rate and transfer equation. Chapter III 
presented a solution of a two level maser of spherical geometry which 
couples to a multilevel set of rate equations; this chapter will 
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describe how those equations can be used graphically to understand 
the characteristics of the maser action and to predict the intensities 
of proposed maser emissions. 


B. The Parameters of the Maser Emission 

The formulae of Chapter III relatesthe pumping rate; the saturation 
state^ and the intensity to the physical conditions of the maser region. 
The parameters are used in a dimensionless form and are scaled with 
the physical conditions in the following manner. The physical diameter 
is scaled to a dimensionless radius by: 


* 


R = C D(l_r) 


where D is the physical size in centimeters; C is the scale length: 




( IV-B-l) 


and r is the unsaturated population ratio n^/n^; so that (l-r) is 
the population difference An/n^. If there is an unsaturated .core in 
the region, its size is similarly scaled by: 

a = C L(l-r) 


where L is' the physical size of the core in centimeters. The pumping 
rate is given by (l + T/A) where A is the spontaneous emission 
coefficient. The background and emission intensities are given by: 


3 




) 


T 

o 


T B = 


Ic 

2kv £ 


I c 
o 


2kv c 


( IV-B-2) 
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where is the emission brightness temperature and the background 

brightness. The state of saturation of the region is given for partially- 
saturated clouds by the size of the unsaturated core; the boundary 
between the core and the saturated shell is defined as the point where 
the spontaneous and the stimulated transition rates are equal. If there 
is no core and the stimulated transition rate exceeds the spontaneous 
rate throughout the cloud, the emission is characterized by the satura- 
tion ratio K of the stimulated transition rate divided by the spontaneous 
rate at the center of the cloud. 

C. A Graphical Sol ution for the Emission 

Equations . II I-B-10, III-C-15, and III-D-3 give the peak brightness 
temperature for a spherical cloud, but are hard to use except for the 
case of unsaturated emission because of the interdependence of the 
intensity and the saturation state. This section will present a two 
step graphical analysis of the emission which additionally provides a 

4 

clear picutre of the character of the maser amplification and the sensi- 
tivity of various parameters. The state of saturation in the cloud will 
be considered first, followed by the intensity determination. 

For a region with an unsaturated core the size of the core is re- 
lated to the total size and the pumping rate by (lII-C-l4); 
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The region is completely saturated when a = R in the above equation, 
and has just reached complete saturation when a = 0. For regions with 
a higher saturation level the appropriate equation is (lII-D-2); 



(IV-C-2) 


If the scaled size CD 'is chosen as the independent variable, 
the pumping rate becomes the dependent function, and the population 
and saturation ratios become parameters of the relationship. Equations 
(iV-C-l) and ( IV-C-2) then divide the size versus pumping rate field 
into regions of unsaturated, partially saturated, and over saturated 
amplification, with the boundaries dependent on the parameters. Figure 
4 gives the boundary between unsaturated and partially saturated emission, 
and figure 5 ‘the boundary between partial saturation and over saturation, 
with zero background intensity and population ratios between 1 and 10 . 

In figure 4 the unsaturated region is to the left of the function line, 
and in figure 5 the over saturation region is to the right. In the 

ft 

region between the two curves, for a fixed population ratio, the cloud 
has 1 an unsaturated core. For a population fraction less than 0,01 the 
partially saturated region at small radius values is negligible, so that 
the emission can change rapidly from unsaturated to fully saturated 
with small changes 'in size or pumping rate. 

If the pumping rate and the size of the cloud meet in a region of 
over saturation, the saturation ratio is an important factor in determining 
the emission intensity. Figures 6 and 6a show the lines of constant 
saturation for two values of the population ratio. In both figures, the 


unsaturated limit is also shown to indicate the limited conditions which 
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can produce partially saturated emission. For values of the saturation 

3 

parameter greater than 10 the boundary moves very little because the 
saturation is depleting the upper state population so strongly that 
further increases in the intensity can have little effect. (Recall 
that the population ratio r is the value produced by the pumping 
process without stimulated emission, not the local value of the popula- 
tion which varies radially in the cloud.) In figures 6 and 6 a the value 
of 0.01 for (l-r) corresponds to a one percent population inversion, 
and a value of 0.1 to a size percent inversion. This is the range of 
inversion needed to explain the microwave emission. 

The effects of the background intensity on the saturation limits 
can be seen in figure 7 for a population inversion ratio of 0.01, The 
background radiation increases the degree of saturation for small 
pimping rates and dominates the saturation level until the' size of the 
region is sufficient for the spontaneous and stimulated transitions in 
the region to affect the emission. For microwave emission the line 
temperature T^ is 0 . 08 °K so even the isotropic 2 . 7 °K cosmological 
background controls the saturation state at small scaled sizes. For 
infrared lines around 100 microns the line temperature is l43°K; at this 
wavelength only a rather bright background would affect the saturation 
state. The appropriate background temperature to use is the brightness 
temperature of the background emission multiplied by the fraction of the 
hemisphere the background source fills as seen by the rear of the maser 
cloud. 

Equation (iV-C-l) also gives a relationship between the overall size 
of the region and the size of the unsaturated core. For partially 
saturated regions the core size determines the degree of saturation of 
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the maser emission. Figures 8 and 9 show the relationship between the 

core size and the total size for zero background and population ratios 

of 0.01 and 0.10 respectively. The difference between the unsaturated 

and the saturated core sizes for a fixed pumping rate ,is less in figure 9 

because the greater inversion allows the intensity to increase as a 

faster rate with increasing size, hence bringing on saturation at a lower 

value of C D. The most interesting characteristic of these two figures 

is the rapidity with which the core disappears after a critical size is 

reached. In figures 8 and 9 the critical size occurs at the elbow of 

the curves, with the core size going to zero with only a small change 

in the total size. The critical size of the core is seen to vary only 

slowly with increasing values of the pumping parameter. If the OH or 

the H^O emission regions are only partially saturated with a core near 

this critical limit and the region is subject to turbulence or other 

* 

variations which could change C D the saturation state and the intensity 

could vary greatly over short periods of time as the core emission of 

exponential growth flickers on and off. For a background intensity with 

T q /T less than 1000 the relationship between core size and total size 

is nearly unchanged; for higher values the curves of figures 8 and 9 

* 

move slightly to the left, with the critical value of C D about fifty 

percent lower for T /T = 1000, 

o v 

Once the saturation state and the parameters of the maser region are 
known, the determination of the emission intensity is straightforward. 
Equation (lII-B-lO) gives the maximum possible brightness temperature, 
from an unsaturated region as: ' 
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(iy-c-3) 



For the paz-tially saturated case the intensity is (lII-C- 15 ); 



If the region is over saturated,, the saturation parameter k characterizes 
the emission and the intensity is given by (lII-D-3): 




1 ' 

-1 . 

( 1 V-C- 5 ) 


The intensity of the spherical region and its dependence on the 
saturation state is presented in figures 10 and 10 a for population inver- 
sions of one and six percent. The unsaturated limit is the upper limit 
for the emission from a maser region of the given parameters, and the 
inversion limit is the lowest emission which can" occur ■with any popula- 
tion inversion. For scaled sizes over a few hundred the saturation state 
is the most critical parameter of the emission intensity, with more 

than 10 orders of magnitude between the unsaturated intensity and the 

* 

inversion limit at C D = 1000. As the population ratio decreases the 
maximum intensity approaches the inversion limit. In figure 11 the 

unsaturated maximum intensity is shown for population inversions down to 

-3 

10 percent. This can be used to set upper limits to the possible 
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emission intensity even if the saturation state is not known. The lower 
limit on the intensity for any population inversion can be seen in 
figures 10 or 10a and is independent of the population ratio. 

t 

D . Summary 

Figures 4 throxigh 11 give a straightforward method for determining 
the emission intensity from a spherical region of population inversion. 

The pumping model gives the population ratio and the pumping rate, and 
the physical conditions gives the scaled size. Figures 4 to 9 are used 
to determine the state of saturation of the region, and figures 10 and 
11 then directly give the peak brightness temperature. The average 
brightness temperature can then be determined from (lII-E-1, 2 or 3) and 
the intensity is given in watts pel' square centimeter per hertz per 
sterradian by: 

2k ( T u) / P 1 l \ 

I = (watts cni Hz*str'* J ' y (iV-D-l) 

X 

Finally the flux can be found from: 

8vk /t \ /Cl \ , 0 

? v = g ( " S 4rt ~ " / ( watts y (IV-D-2) 

X 

Then the collecting area and the receiver bandwidth are all that is 
•needed to determine the pouter received at the detector. 


67 



Figure 
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Figure 


CAPTIONS FOR FIGURES 4 THROUGH II 

4. The boundary between the unsaturated state and the partially 
saturated state for spherical masers as a function of the 
pumping rate and the scaled size, for selected parameters 

of the population fraction. The size is scaled by 
equation (iV-B-l), 

5. The boundary between the partially saturated state and the 
fully saturated state for spherical masers 'as a function 
of the pumping rate, the scaled size, and the population 
fraction. 

6. The saturation limits of a spherical cloud with a population 
fraction of 0.01, showing the region of partial saturation 
between the unsaturated limit and the saturation limit, and 
indicating the degree of oversaturation with lines of 
constant K. 

6a. The saturation limits for a population fraction of 0.10„ 

7. The effects of the background intensity on the partial 
saturation boundaries as a function of the ratio of background 
brightness temperature to the line temperature hv/k for 

a population fraction of 0.01. 

8. The relationship between the size of the unsaturated core 
and the total size for zero background intensity and a 
population fraction of 0.01, at various pumping rates, showing 
that the unsaturated core disappears with a limited variation 
in the total size. 

9. The relationship of figure 8 for a population fraction of -0.10. 

10. The maximum brightness temperature at the center of the 
observed disc for a spherical maser with a population fraction 
of 0.01, as determined by the state of saturation of the cloud. 
The unsaturated limit is the highest possible brightness for 
this population fraction, and the inversion limit is the 
minimum; the saturation limit is the intensity which just causes 
the unsaturated core to disappear. 
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Figure 10a. The central brightness of a spherical cloud with a popula- 


tion fraction of 0.10. 

Figure 11. The maximum central brightness of spherical masers for a 
range of population values as a function of the scaled 
size of the region. 
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Figure 4, Unsaturated State Limit 
70 
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Figure 5. Complete Saturation State' Limit 
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Figure 6. Oversatir ation State Limits 
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Figure 6a. Over saturation State Limits 
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Figure 10. Peak Intensity for Population Fraction 
( 1-r) = 0.01' ‘ 
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Figure 10a. Peak Intensity for Population Fraction 
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Figure 11, Unsaturated State Peak Intensity 
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CHAPTER V 


POPULATION OF THE OH ENERGY LEVELS 


A. Introduction 

In the decade since the discovery of the OH emission a number of 
mechanisms have been put forward to explain the maser pumping. The 
four mechanisms that will be discussed here cover the major exciting 
pathways available in the OH molecule, though the actual pumping which 
occurs in astrophysical masers may be a complicated combination of any 
or all of these mechanisms. The first OH emission discovered was 
associated with the ionized gas near hot stars. This association suggested 
that ultraviolet radiation could be the exciting power for the masers 
through resonant lines near 3080 A. Litvak et al, (l 966 ) and Perkins 
et al. (1966) proposed a UV absorption mechanism which required selective 
absorption in the UV lines to provide the microwave inversion. Infrared 
astronomers have also associated some OH sources with infrared stars, 
which can pump the OH through the ' 2 . 8 and 1.4 micron vibrational bands, 
or for lower temperature infrared clouds the pumping could be through the 
100 micron rotational levels. Studies of these pathways have been 
presented by Litvak (1968, Litvak et al. (1969)* anc * a £ a i- n by Litvak 
( 1973 ) and by Shklovsky (1967). Townes ( 197 ^) and others have also 
proposed that the OH emission could be pumped by the collisional action 
of atoms or molecules with the OH. Most 'of the other mechanisms proposed 
can be lumped under the title of chemical pumping where the OH is created 
in an excited state by the formation process. This later is so dependent 
on the assumptions of the formation mechanism that any desired popula- 
tion distribution could be produced with careful assumptions, and is 
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hence too vague to provide further insight. We have made numerical 
calculations of the populations produced by the other mechanisms and 
will present them here. Our desire is to examine the far infrared 
emission which would result from the pumping, and consider the micro- 
wave emission only as it can help -determine the physical conditions of 
the maser region. This allows us to ignore the hyperfine structure of 
the OH molecule without significant loss. The energy structure of the 
ground vibrational state is shown in figure 20 for reference, with the 
dipole allowed transitions and their wavelengths in microns. 

In the derivations of Chapters II and III the interactions between 
the transfer equation and the population equations have all been folded 
into a modified transfer equation, whose solution was then presented 
in Chapters III and IV., The population solution of the emission problem 
then becomes an eigenvalue solution of the matrix of transfer coefficients 
and can be written as the system of equations: 


dn. 

g i x dT 


E W..n. - n. V'W. . 
Jl J IjLu 13 

j J 


= 0 


where 


(V-A-l) 


W. . = A. . + B. .1 + C. . 
xj xj ij 13 


and the n.^ are the normalized populations corrected for degeneracy. 
The solution of the, above matrix will then give the population fraction 
as n^/g^ since this is the value of the population which will be 
needed in the transfer equation. 
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B. Ultraviolet Pumping 


The ultraviolet pumping mechanism proposed by Perkins et al. ( 19 66 ) 

a 

and refined by Litvak et al. ( 1966 ) involves the absorption of 3080 A 
radiation through the transitions from the ground state lambda doublet 
of the OH molecule to the six lowest rotational levels of the first 
excited electronic state of OH, followed by spontaneous decay to the 
various rotational and vibrational levels of the V = 0, 1, and 2 
ground electronic state. The exciting source is’ assumed to be an 05 
type star at a distance of 2 light years. The only radiation absorbed 
is through the six transitions leaving the lowest Lambda doublet 
although transitions involving the other rotational levels of the 
ground electronic state are nearly the same wavelength. This selective 
action could occur by resonance with' some line structure in the 05 star, 
or through the filtering action of an intervening gas such as hydrogen. 

Figure 12 shows the populations of the OH states which results 
from this mechanism for various distances of the exciting source. The 

■ I 

ground state lambda doublet is not inverted, so this simple mechanism 
cannot provide the observed microwave emission. However, the lambda 
doublets of two excited *11 ^ states are inverted and could lead to 

detectable maser emission if the exciting source were closer than two 
light years. If the exciting' source were very close, the 12-10, 11-9* 
8-4, and 7-3 inversions could lead to observable infrared emission at 
wavelengths near 100 microns. If the exciting source were farther away 
the OH cloud would have to be extremely large to give observable 
emission. Unless the exciting source were very restricted in its inten- 

e 

sity outside the 3080' A absorption band the OH would probably be 
disassociated by the star at distances closer than a few tenths of a 


light year. 
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As the ultraviolet radiation propagates through the OH cloud, 
absorption will alter the intensity distribution in the six UV lines. 
Litvak et al. (1966) showed that this selective absorption leads to 
a limited region of population inversion in the observed ground state 
lambda doublet. Figure 13 shows this effect for the UV mechanism 
with the 05 star at a distance of two light years, with the optical 
depth scaled to the strongest UV absorption line. At this distance 
the higher rotational levels have small population fractions, but the 
microwave inversion of the excited lambda doublets is maintained or 
increased over a greater distance than the ground state doublet. Despite 
the low populations it may be possible to detect this emission if the 
region is large enough. This microwave inversion would also not disappear 

i * 

if the exciting source were much closer, and in that case there could be 
detectable infrared emission.' . 

C. Collisional Pumping 

In equation (v-A-l) a term Ch ^ has been included for the colli- 
sional transfer rate. Dipole collisions between OH and the common 
charged particles of the interstellar gas have been studied by Barrett 
and Lilley ( 1957 ) and further by Rogers and Barrett (1968) and Goss 
and Field (1968). However, the regions of OH emission are replete with 
complicated molecules, and charged particles can be uncommon members of - 
the OH medium. 

Townes (197O) and Goldsmith ( 197 ^) suggested that molecular colli- 
sions occurring in the interstellar medium could lead to excitation of 
maser transitions in ammonia and carbon monoxide, and possibly in OH. 
Townes and Gwinn (see Townes 197O) suggested that if the colliding 



particle is a diatomic hydrogen molecule the OH could be excited in a 
preferred rotational state which would lead to the microwave inversion. 
Goldsmith found that in the case of CO the collisions with charged 
particles, although more frequent than molecular collisions, would 
produce only a thermal population distribution because the collisionally 
induced transitions were restricted to dipole allowed changes. Even 
the quadropole transitions produced a nearly thermal population distribu- 
tion. However, as Goss and Field (1968) and others have' shown, there 
are no selection rules for molecule-molecule collisions because the 
interaction forces are very short range. Goldsmith then applied constant 
deexcitation cross sections with exciting cross sections given by: 


-AE 

g.cr. . = g.cr.. exp — - 
1 lj 1 p RT 


j > i 


( V-C-l) 


to find population inversions in the CO molecule. Figures l 4 and 15 

show the effects of this pimping scheme ‘on the OH energy structure, 

-15 2 

using a constant de-exciting cross section of 10 cm , with the 
reverse exciting collisions using the above cross section. There are no 

Q 

selection rules in this scheme, and for densities below 10 Hg molecules 

3 

per cm the spontaneous transition rates of the infrared lines exceeds 
the collisional rates. Hie action of this mechanism is then to excite 
the OH to all possible levels, the variation in rates among the levels 
being only the Boltzmann factor in ( V-C-l), with the excited molecules 
returning to the lower states by spontaneous emission. Since the exciting 
and the de-exciting rates are not symmetrical the OH population distri- 
bution can take on anomalous inversions just as it would if the OH were 
formed in the excited states. For densities below a few particles 
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per cubic centimeter the collision rates are not sufficient to excite 
the OH molecule. As the kinetic temperature increases the exciting 
and de-exciting collision rates becomes more nearly equal, and the 
population distribution takes on a thermal character. For the kinetic 
temperature of 100°K only the l8 centimeter microwave line is inverted 
sufficiently to produce detectable emission unless the Hg density 

g 

exceeds 10 particles per cc. In figures l4 and 15 the populations 

of the odd numbered states except 1 are identical to the populations of 

their even numbered doublet partner, so only the even numbers are given. 

7 -3 

Above 10' cm densities the infrared lines from levels 12 to 10 and 

8 to 4 are sufficiently inverted to produce maser radiation. At a 

kinetic temperature of 1200°K the collisional rate across the lambda 

doublet is nearly equal- in both direction, and the microwave inversion 

6 8 “3 

is disappearing. Between 10 and 10° cm” the inversion of the 

micron line between levels 8 and 4 disappears. The 12 to 10 inversion 

5 -3 

is significant at densities above 10 cm . The higher vibrational and 
electronic states have been included in this model, but the energy 
differences are so large that the collisional pumping from the v = 0 
levels is negligible, and the populations are near zero. 

D. UV and Collisional Pumping Combined 

Figure 1 6 shows the combined effects of the UV pumping scheme 
(without absorption) with the exciting star at two light years, with the 
molecular collision scheme. Because the UV pumping at this distance is 
so weak, the collision effects dominate the excited state populations. 

In fact, the collisions can still produce the microwave inversion if 

3 -3 

the density of Hg particles exceeds 10 cm . Because of the weakness 
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of the UV pumping the infrared emission is the same as in the 
collisional case. 

E. Near infrared Pumping 

The vibrational bands of OH at 1.4 and 2.8 microns can be excited 
and produce maser pumping by resonance absorption of infrared -radia- 
tion. A number of the other molecules in the interstellar medium have 

emission in the range of these transitions, and for temperatures of 

• / 

a few thousand degrees the black body peak occurs between the two 

wavelengths. If the OH is near a star or protostar of this temperature, 

it would be possible for the vibrational bands to be excited. If 

there is a significant amount of dust in the vicinity, the longer far 

infrared portion of -the spectrum could be strongly absorbed, so that 

only the vibrational transitions are excited. If the radiation incident 

on the OH cloud is of equal intensity in the 1.4 and 2.8 micron bands 

« 

and nearly zero at other wavelengths the population will have the 
distribution shown in figure 17 . ' The temperature of a black body with 
equal intensities at 1.4 and 2.8 microns is 2650 °!!, so this is the listed 

* f 

exciting temperature, although the excitation is more correctly a pair 
of equal delta functions at 1,4 and 2.8 microns with a line width equal 
to the vibrational bandwidth. The intensity of the pumping radiation 
is the black body intensity at 26 c yO°K multiplied by the fraction of the 
sky th'at the OH cloud sees filled with the pump source. This is the 
equivalent to the distance function of the UV pumping case, but here 
it is difficult to estimate the physical size of the infrared source. 

For vibrational excitation the parity selection is fixed for all the 
transitions, and there is no parity assymmetry which can produce a lambda 
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doublet inversion. Litvak and others have proposed that the trapping 
of the far infrared radiation by resonant scattering in the OH region 
as the mechanism for producing the microwave inversion. Selective ampli- 
fication of the magnetic hyperfine structure caused by the presence of 
a magnetic field, or a plasma propagation effect such as parametric 
down conversion discussed by Litvak ( 197 O) are needed to produce the 
microwave inversion with the correct polarization properties. 

The near infrared pumping mechanism does produce some very 
interesting results for the far infrared emission pattern. Figure 17 
shows the population solution, with only the odd numbered levels shown, 

I 

as the even numbered levels have identical populations to their doublet 

counterparts; the dashed lines are the populations of the rotational 

levels in the first vibrational state. If the pump source is a large 

part of the sky at the OH region, the lowest ^ eve ^ ^ as a ^ ar 6 er 

population than the ground state, and a very strong infrared maser emission 

4 

at 82 microns is expected. A sky fraction below 10 for the pump source 
does not produce any infrared inversion, and 'as the intensity increases 
the infrared inversions appear and grow more intense. When the pumping 
source fills half of the sky at the OH cloud all the infrared transitions 
connecting the *^ 3/2 an< * ^ 1/2 s ^ a ^ es have inverted populations. 

More intense pumping can lead to population inversions in some of the 
2.8 micron transitions connecting rotational levels of the v = 1 
and the v = 0 states. 
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F. Far Infrared Pumping 


The assumption in case E that the far infrared radiation from the 
pumping source is negligible is not easily justified. Unless the ex- 
citing radiation is a resonant line from some other emission source, 
the radiation in the 100 micron region of low temperature black bodies 
can be important in the OH pumping. Litvak and others have proposed 
that the trapping of far infrared radiation in the OH cloud or from a 
nearby protostar is important in producing the microwave emission from 
some of the Class II OH sources. This section deals with a simple far 
infrared model to show that for some conditions the FIR radiation is 
indeed important, at least to the production of infrared masers. 

■•The pumping proposed here involves the absorption of black body 

radiation from low temperature objects in the vicinity of the OH cloud 
K 

as a function of the sky fraction they cover. All the infi’ared transi- 
tions of the three lowest vibrational states and the vibrational transi- 

t 

tions among these three states. For temperatures of a few hundred degrees 
the black body peak, and hence the greatest transition rate is near the 
100 micron wavelengths of the rotational transitions. The population 
distribution among the low rotational states should be nearly thermal 
with little infrared inversion possible. The only source of inversion 
pumping for these very low temperatures is the dilution of the pumping 
radiation by the background space. If the pump only occupies a small 
fraction of the sky however the intensity will be too weak to provide 
detectable emission lines. Since the absorption of black body radiation ■ 
has no parity selection rules there will be no microwave inversions, 
with the populations of the doublet pairs identical. If the temperature 
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of the black body pump is in the range of a few thousand degrees, the 
vibrational transitions at 1.4 and 2.8 microns are strongly excited, as 
in the near infrared pumping model of Section E. 

Figure l8 shows the population distribution which results from the 
absorption of 500°K black body radiation. In order to see the infrared 
emission the OH cloud must be nearly enclosed by the pumping source, and 
then only weak inversions are produced in the 5 to 1 and higher transi- 
tions. In order to be detectable the infrared emission resulting from 
this pump will have to come from a very large OH cloud. Figure 19 shows 
the population distribution produced by a 2000°K pumping source. Pumping 
to the vibrational levels and the ensuring decay to the excited rotational 
levels affects the ground state population even for small sky fractions 
allowing the 5 to 1 trahsition to become inverted at a sky fi'action of 
0.1. Sky fractions above one percent produce 'weak inversions in the 
11-9 an d 15-13 lines, and larger pumping rates produces infrared inversion 
in all the AJ = +l transitions, ' and the 50 micron band of AJ = 0 
transitions are nearly inverted when the sky in ninety percent filled with 
the pumping source. For high pumping rates the populations of the v == 1 
states begins to become significant. The population distribution of 

* * v 

figure 19 is not greatly different from the near infrared pumping model 

-2 

of figure 17 for sky fractions above 10 if the difference in temperature 
is considered. 

G. Summary 

The four pumping mechanisms discussed in this chapter cover the 
most likely exciting pathways available to the OH molecule which do not 
require the creation of OH in an excited state. Under two of them the 
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microwave main line inversion is easily created, and under all of them 
there are some infrared line inversions. A more detailed analysis of 
the energy and transfer structure of the OH molecule which includes 
the hyperfine structure is necessary to model the class II satellite 
line emission, though efforts have been made in this direction in the 
papers mentioned by Litvak and others. The ease with which the colli - 
sional pumping model produced the microwave inversion implies that 
collisions can indeed be important in the emission of the interstellar 
molecules, a view that has long been held by Townes and others. The 
common occurrence of infrared inversions means that infrared observations 
can make significant contributions to the study of the interstellar 
molecules . Current observing capabilities in the infrared region around 
100 microns are limited in sensitivity, and the next chapter will consider 
the possibility of detecting the infrared maser emission of the various 
pumping models . 
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Captions for Figures 12 through 20 


Figure 


Figure 


Figure 


Figure 


Figure 


Figure 


Figure 


12. The population fraction resulting from the UV pumping 
model including vibrational levels as a function of the 
intensity or distance of an 05 exciting star. 

13. The population distribution of the UV pumping model including 
* the effects of attenuation of the UV radiation by the OH, 

as a function of the optical depth in the P^l) transi- 
tion at 3080 A. 

14. The population distribution of the collisional pumping model 
- with the kinetic temperature of the H^ col-liding particles 

at 100°K. The cross section of the collisions is assumed to 

-15 2 

be 10 J cm . 

15 . The population distribution of the collisional pumping model 
for a kinetic temperature of 1200°K, which corresponds to 
the thermal line width of the microwave emission. 

16. The population of the OH molecule with the combined action 
of UV pumping by an 05 star at a distance of 2 light years 
and collisional pumping at a kinetic temperature of 100°K, 

17 . The population effects of near infrared pumping, where the 
exciting radiation is confined to the 1.4 and 2.8 micron 
bands, with the ratio between the two band strengths corres- 
ponding to a temperature of 26^>0°K. The dashed lines are 

the populations of the first vibrational state. The intensity 
of -the pumping radiation is given by the black body intensity 
at 265 O °K times the fraction of the sky that the pump fills as 
seen by the OH cloud. 

1 8 . The population distribution of the OH molecule with pumping 
by a 500°K black body filling a fraction of the sky. Only 
the odd numbered levels are shown, with their even numbered 
counterparts having the same population. 
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Figure 19 .- The population distribution for black body pumping at 
2000 °K as a function of the dilution factor. 

Figure 20. The energy level structure of the lowest vibrational state 
of the OH molecule showing the allowed transitions and 
their wavelengths, and the numbering scheme used in 
figures 12 through 19 . 
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CHAPTER VI . 


THE INFRARED INTENSITIES OF OH MASERS 

A. Introduction 

The discovery of the OH emission in 1 965 stimulated a steady stream 
of theoretical discussions which have attempted to explain some or all 
of the important characteristics of the OH maser emission. The problem 
of polarization of the OH emission has yielded least to the efforts in 
part because the polarization is a characteristic of the propagation of 
the emission as it passes through the cloud and is dependent on the 
pumping mechanism of the OH. Previous chapters of this work have 
analyzed several of the important pumping mechanisms available to the 
OH in interstellar space and along with the papers published by other 
authors have shown that the observed OH intensities can be explained by 
one or another of these pumps. To make further headway against the 
problems of polarization, and observed size, and of the conditions in the 
OH regions, a determination of the active pumping mechanism for each 
emission source must be made. The results of the pumping analysis of 
Chapter V indicates that population inversions can occur in the infrared 
transitions of the OH molecule in different proportions depending on the 
active pumping mechanism involved. This implies that successful obser- 
vation of the infrared line emission from the Oil sources can aid in 
determining the pumping mechanism and physical conditions of the region. 

Using the model of the maser emission in spherical regions presented 
in Chapters III and IV we can describe the conditions of the emitting 
region as derived -from the observed microwave intensities under various 
pumping mechanisms, and apply these conditions to the prediction of the 
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infrared line emission. As shown in Chapter V a pumping mechanism which 


allows parity effects., such as the TJV pump or one of the collisional 
models, and which produces a ground state population inversion also 
produces population inversion for some infrared lines. Of course mecha- 
nisms which require propagation effects to produce the microwave inver- 
sion can also have infrared inversion, but we cannot use the microwave 
observations to analyze the conditions of the OH regions in these cases 
using our simple transfer model. For these pumping mechanisms a set of 
physical conditions consistent with the observed emission can be assumed 
in order to estimate the resulting infrared emission. Section D will then 
consider the detectability limit of current and projected infrared 
instruments and the limits this sets on the possible OH emission. 

B. Characteristics of OH Emitters 

The complicated nature of the OH emission has been discussed in 
Chapter Ij here we will mention the properties of the emission in the 
sources W3 and Orion OH which are appropriate . to the models of the 
maser amplification and pumping mechanism discussed in Chapters III and 
IV. 

The W3 region of OH emission is one of the -most heavily studied 
and complicated of the microwave OH emitters. W3 is in a HII region, 
which is maintained by one or more stars of spectral class 0. Infrared 
observations have also revealed significant near and far infrared emission 
similar to hot dust 'and IR stars in the vicinity of the OH microwave 
sources, and a number of other complex' molecules have been detected in 
the region. The physical sizes observed may only be vaguely related 
to the actual dimensions of the OH emitters, but they indicate a number 
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of small regions with diameters between 10 and a few hundred astrono- 
mical units (A.U.) scattered over a region of about 2000 A.U. The 
lines observed in emission are at frequencies of 1 665 and 1667 MHz, with 
line widths of 5 to 15 kHz. The effective brightness temperatures, 
given by: 



*5 o 

have steadily increased from 10 K to the presently accepted values 
11 13 

between 10 and 10 degrees as the resolution of long baseline inter - 

, -8 _g 

ferometry has lowered the observed size to 10 to 10 sterradians. 

Other observations of OH absorption clouds indicated that the OH eoncen- 

3 -3 

tration is between 1 and 10 cm and the densities of other molecules 

IQ -3 

and atoms less than 10 cm . 


The OH source in Orion is located near the trapezium stars, an 

1 

area generally associated with star formation. As in W 3 there is signi- 
ficant observed infrared emission from a large region surrounding the 
OH source and an infrared star of about 2000 °K nearby. The optical 
obscuration in the region has been estimated at 20 to 100 magnitudes 
due primarily to a large amount of dust and cool. gas. The physical size 

observed is less than 2 X 10 cm, and the linewidth at 1665 and 1667 MHz 

* 

is about 2 kHz. 


C. The Infrared Emission of the Orion and W 3 OH Sources 

For the ultraviolet and collisional pumping models discussed in 
Chapter V there is significant inversion of the microwave energy levels 
which would lead to the maser amplification of the 1 665 and 1667 
MHz lines. In these models we will use the .observational properties 
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mentioned above to estimate the physical characteristics of the OH clouds, 
and then use these estimates to calculate the expected infrared emission. 

Hie other pumping models do not show inversion in the lambda doublet 
which gives rise to the 1 966 and 1667 MHz lines, but propagation effects 
such as resonant scattering infrared trapping, or parametric conversion 
could lead to the inversion of the satellite lines at l6l2 and 1720 MHz, 
or possibly in the main lines. For pumping models which do not have 
intrinsic microwave inversion we will assume some characteristics consis- 
tent with the microwave observations. 

In order to determine the intensity of a maser -amplified line using 

the relations of Chapter III or the graphs of Chapter IV the following 
* > 

procedure will be used. Chapter II gives the spontaneous coefficient and 
the wavelength of the transition; these are combined with the physical size 
and the OH density to determine the scaled size from (iV-B-l), also using 
the population of the upper state of the maser emission. The numerical 
pumping model calculation which determined the population also gives the 
population ratio and the pumping rate P of the transition using (lI-E-5)j 

3 . 

(II-E-8), and (iI-E-ll), The scaled size, population ratio, and the pumping 
rate are then used with figures 4, 5 , or 6 to determine the state and degree 
of saturation in the spherical region. For partially saturated regions 
figures 8 and 9 permit the determination of the core size using the calcu- 
lated parameters. Finally, figures 10 and 11 determine the scaled intensity. 
The peak brightness temperature so determined is averaged with (lII-E-l), 
(lII-E-2), or (lII-E-3)^ and the intensity can be found with the line 
width and source size using (iV-D-l). 

If the microwave intensity is to be used to determine the parameters 

of the 1 OH cloud, the above procedure is reversed. The pumping model and 

* 

the observations are' used to determine the scaled intensity, which sets 
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range on the scaled size that depends on the saturation state. The 
pumping rate from the pump model is combined with, the population ratio 
and this range of sizes to choose a consistent set of values for the 
size and saturation. Then the density and molecular characteristics 
are used to determine the physical size. The procedure can be iterated 
to verify the chosen values. The parameters so found are then used in 
the foregoing manner to find' the infrared intensities. 


In the transition matrix of the OH molecule the rates connecting 
the J levels of the ^2./2 s l a l es are greater by factors of 10 or 
more than the rates connecting the s ^ a ^ es ^ 3/2 s ^ a ^ es * 

This leads to population inversions in the crossing transitions under 


most of the pumping models considered^ with the wavelengths and molecular 


constants of: 


wavelength 

A.. 

0 1^4 
11 

b? 

82.6 p 

3.5 x 10 2 sec 1 

174°k 

96.4 p 

8.9 X 10 ^ sec ^ 

i49°k 

115.1 d 

-3 -1 

3.1 X 10 sec 

! 

125 °K • 

134.9 p 

-3 -1 

1.2 X 10 J sec 

• 107°K 

18 cm 

-11 -1 
7.11 X 10 sec 

o.o 8 °k 

6.3 cm 

l.l 6 X 10"9 sec -1 

0 . 22°K 


1. Ultraviolet Pumping 

In order to produce the microwave inversion with ultraviolet 
pumping the selective absorption of the UV radiation propagating into 
the cloud must be considered. Hie population distribution with this 
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absorption included is not significantly different for the infrared lines 
than when absorption is not included. The microwave inversion occurs 

i 

over a limited distance, while the infrared inversions are maintained 
over significantly larger distances. With the infrared populations 
which occur under pumping from 2 light years the increased size of their 
inversion regions is not of importance, since the emission will remain 
below detectable levels. If the exciting source is moved closer, to 
1000 A.U., the populations of the infrared transitions are much larger, 
and the population of the 6.3 cm .transition has increased to detectable 
levels . Of course at this close distance the OH dissociation rate could 
increase to a degree which would destroy all the OH emission. If this 
dissociation is tempered by a high formation rate, the entire emission 
complex of the W3 region can be explained by the selective absorption 
creating a shell of microwave inversion.. Density variations or turbulence 
in the OH cloud surrounding the exciting star could then lead to. the 

t 

formation of numerous spots of microwave emission, as is observed. However, 
the inversion in the 6.3 cm line of the ^2./2 f ^ ~ doublet could 

occur over the entire region. Since the population in this line is an 
order of magnitude below the l 8 cm populations and the pumping rate in 
this level is lower, the 6.3 cm emission could be lower in intensity, as, 
is observed. The large extent of the population inversion for the infrared 
transitions leads to significant emission intensities. 

For UV pumping from 1000 A.U. the inversions, their populations, their 
population ratios, and their pumping rates [found by numerical solution 
of equations (lI-E~l), ( 1 I-E- 5 ), (lI-E- 8 ), and (iI-E— ll)], are 
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\ 

V n <H 

(1-r) 

r /a 

a 

6. 3 cm 

4 .o x io ~ 2 

0.32 

3 x 10 9 

96.4 p 

-2 

1.0 X 10 

0.8 

200 

96.5 ft 

O 

6.7 X 10 

0.6 

105 

115.1 p 

- 3 - 

1.02 X 10 J 

0.9 

290 

115.2 p 

8.8 X 10 ”^ 

0.95 ! 

300 


In the UV pumping scheme there are three sources to model; the Orion OH 
source, the \V3 source, and the extended W3 source with the emission 
region ten times larger for the infrared lines and for the 6.3 cm line 
than for the l8 cm line. The scaled size, saturation state, and average 
brightness temperature using figures 4 and 5 to determine the state of 
saturation, (and figure 6 for the microwave lines, which are fully 
saturated), and figure 10 and 11 for the intensity. The set of results 
which consistently fits the l8 cm observations of the Orion OH source are 



* 

C D 

* 

State 


p 

I (w/cm Hz str) 
V 

6. 3 cm 

380 

FS, K = 50 

10 11 

-l 4 

7.0 X 10 

96.4 pm 

8.0 

PS, a/R =0.2 

2 X 10 5 

_i 4 

5.9 x 10 

96.5 pm 

7.2 


4 

9 x 10 

• 

_i 4 

2.7 X 10 


f 


F.S. = fully saturated, the extent indicated by the parameter K. 
P.S. = partially saturated. 

U.S. = unsaturated. 
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The parameters and conditions of the region as determined from the '18 cm 
observations are: 

l4 

MW observed size = 3 X 10 cm 

1 1 }. 

physical size = 8 X 10 cm 
2 -3 

n 0H = 10 Cm 

For the W3 OH source the expected results are: 


X 

* 

C D 

State 

( t b>° k 

2 

I (w/cm Hz st: 
V 

6,3 cm 

500 

F.S. K = 9 

CVJ 

r-i 

O 

H 

-13 

7.0 X 10 J 

96.4 pm 

19.0 

P.S. a/R = 0.2 

g 

1.5 X 10 

4.4 x 10“ 13 

96,5 pm 

11.6 

P.S. = 0.2 ! 

7 x io^ 

-13 

2.1 x 10 J 

115.1 pm 

3.4 

P.S. = 0.8 

4 

3 x 10 

6.3 x i0“ 15 

115.2 pm 

2.3 • 

P.S, ,= 0.8 

’ 4 

1 X 10 

-15 

2.1 X 10 y 


where the l8 cm observations gave 


MW observed size = 9 X 10 

18 

physical size = 1.8 X 10 ^ 


l4 


2 -3 

n 0H = 10 cm 


cm 

cm 


For the extended W3 source the results which match the l8 cm and the ' 
6.3 cm observations (the latter being used to estimate the size of the 
entire OH region while the former determines the smaller l8 cm inver T 
sion region) are: 
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\ in pm ' 

i *. 
C D 

- State 

< t b)° k 

2 

I^( w/cm Hz str) 

96.1* 

16.3 

P.S. a/R = .3 

3.2 X io? 

-12 

9.5 x 10 

96.5 

10.1 

P.S. = .3 

1.0 X 10? 

3.0 x 10 -12 

115.1 

9-7 

P.S. ‘ = .9 

k x io5 

8.k x I0~ lli ' 

115.2 

9-3 

P.S. = .9 

3 x 10? 

6.3 x io~ l4 


with the 6.3 cm and l 8 cm observations giving: 

"oh ” 10 c "' 3 ' - 

15 

physical size = 7-3 X 10 cm 
MW observed size = 1.0 X 10 ^ cm 

2. Colli si onal Pumping 

The pumping of the- OH molecule by low temperature kinetic collisions 
with hydrogen molecules is not a likely model -for the OH emission in the 
Orion vicinity because of the large amounts of infrared radiation with • 
a spectrum similar to a 700°K black body. Further, if the observed line 
widths are due to thermal broadening, then the implied thermal tempera- 
tures are rather higher than, the 100 degree temperature of the cool 
collision model. Nevertheless, this model is able to produce the inver- 
sion needed for the microwave emission, but the infrared emission level 

is below the detectability limit, unless the -density of colliding particles 

8 -3 

becomes greater than 10 cm . For a density equal to that value, the 
population results are: 


X 

V"0H 

(1-r) 

r /a 

a 

96 n 

2.1*5 X 10 -3 

0.05 

90 

115 

_k 

1.5 X 10 

’ 0.3 

70 
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and the emission determined from figure 10 and 6a is: 


\ 

* 

C D 

State 

< T b> 

I_^( w/cm^Hz str) 

9 6 u 

29 

F.S , K = 1 

4.7 x io 4 °k 

t -l4 

1.4 X 10 


If the collisional temperature is increased to the thermal line widths 
1200°K, the emission in the infrared states becomes more visible. The 
population inversion of the microwave transitions is reduced by the de- 

4 -3 

exciting action of the collisions, but for densities above 10 cm the 

8 

microwave intensity is still achievable, and at 10 the infrared line 
emission is significant. The lines and their population parameters' at 
that density are: 



Applying these results to the determination of the intensities and satura- 
tion states with figures 4 through 11, we find for the Orion OH source 
that: 



where the l8 cm observations give 

"oh - 1q2 an ' 3 

x- i4 

physical size = 6.3 X 10 cm 
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For the W3 emission: 


\ 

* 

C D 

State 

< t b> 

2 

I (w/cm Hz str) 
V 

115 n 

10.8 

P.S, a/R = .6 

7.3 x 10 2 * * * 6 * °K 

1.5 x io" 12 

13b- n 

16.3 

F.S. K = 1 

*5 0 

1.2 X Kr K 

1.8 X 10 


using the l8 cm value of : 

n 0H = 1C)2 Cm_3 

1*5 

physical size = 5 . 1 X 10 cm 

In this pumping scheme the extended W3 source is inappropriate because 
the estimated 6.3 cm inversion is very small. 

3. Near Infrared Pumping 

The unusual near infrared pumping scheme • described in Chapter V 

easily produces infrared population inversions. However, since there 

are no parity effects included in- this model, there is no microwave 

inversion to match the emission with. If we make some assumptions about 

the size and density of the OH sources, the infrared emission can be 

predicted. The assumed diameter of the region is chosen to be three 

times the observed microwave size, to account for reduction in size 

expected in the microwave emission. Since OH absorption studies of the 

W3 region show a /n ~ 10 and other 'work finds n IT 10° cm - the 
UK n H 

2 -3 

value of n n „ = 10 cm will be used for this pumping model. The same 

UK 

parameters will be used for the far infrared pumping model which also does 

not produce a microwave inversion. If the twin delta function pumping 

source has an effective brightness temperature of 2650 °K the infrared 
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_2 

populations reach a maximum for a sky fraction of Q£/2rt = 10 . For 

further increases in the pumping intensity the populations of the 

vibrationally excited state become significant, with little change in 

the v — ‘ 0 populations. The population results for a sky fraction of 
-2 

10 are : 


X 

n / n 0H 

(1-r) 

r /a 

a 

82 p 

7.6 x io -2 

.6 

93 

96 li 

1.8 x 10" 2 

.4 . 

100 

115 n 

8.2 X 10" 3 

.4 . 

560 

135 n 

3.0 X 10~ 3 

.3 

1090 


These lead to the following intensities and conditions for the Orion OH 
source: 


X 

X 

C D 

State 

( T B> 

2 

I^(w/cm Hz str) 

82 p 

6.1 

F.S. K = 1 

2.0 X 1 0 6 °K 

8.1 X io -13 

96 p 

8.9 

P.S. a/R = 0.4 

1.0 X 10 6 °K 

3.0 x io -13 

115 p 

9.1 

P.S. a/R = 0.4 

7.2 X io 5 °K 

-13 

1.5 X 10 J 

134 p 

9.0 

P.S. a/R = 0.5 

1.0 X io 5 °K 

-13 

1.5 X 10 


With the chosen parameters of : 

* i 4 

physical size = 6 X 10 cm 

2 -3 

n 0H = 10 cm 
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For the W 3 source the results are: 


X 

* 

C D 

State 

< T B> 

I^( w/cm^Hz str) 

82 [i 

11.2 

P.S. a/R =: 0.5 

1.5 X 10 7 °K 

6.1 x io~ 12 

96 p 

13 • 

P.S. a/R = 0.1 

9 X 10 6 °K 

-12 

2.7 X 10 

115 H 

11.7 

P.S. a/R = 0.7 

1.4 x 10 6 °K 

-13 

2.9 X 10 0 

13^ p 

10.6 

u.s. 

9.5 X 10 5 °K 

1.4 x io“ 13 


using: 


physical size = 6 X 10 "^ cm 

in 2 '“3 

n 0H = 10 0151 


4 . Far Infrared Pumping 

For the OH source in the Orion region the far infrared pumping 
mechanism can be appropriate for an exciting temperature of 500° K - 
The presence of the HII region in- the vicinity of the W 3 source might 
"indicate that the exciting temperature of 2000°K could be a more valid 
estimate. The infrared emission in the trapezium region of the Orion 
source could also produce 2000°K exciting radiation for the far infrared 
pump. The lower pumping temperature produces a lower total intensity ^ 
so that the populations of the infrared transitions is not sufficient 
unless the pimping source occupies most of the sky as seen by the OH 
cloud. In this instance the population results for the Orion OH source 
are: 


Ilk 




n 

u 

( l-r) 

r /a 

a 

96 p 

1.4 X io~ 2 

-3 

0.06 

76 

115 p 

7.0 x 10 

0.27 

290 

135 H 

3.5 x 10' 3 

0.3 

130 


and the emission from Orion OH is: 


\ 

* 

C D 

State 

< T B> 

p 

I^(w/cm Hz st: 

96 [i 

43 

P.S. a/R = 0.3 

6 X 10 5 °K 

1.8 X 10" 13 

115 p 

23 

F.S. K = 1 

4 O 

3.2 X 10 K 

6.7 x io" 1 ^ 

135 p 

20 

F.S. K = 3 

1.6 X 10 3 °K 

2.4 X 10 "' 16 


where the parameters used are: 


5 i. A f 

X 10 cm 

n OH = 10 cm 


At the higher piunping temperature of 2000°K population inversions in 

the infrared lines become detectable for sky fractions above ten percent. 

For a sky fraction of 50 percent the far infrared lines will be strongly 

inverted, and the pumping rate sufficient to keep the infrared emission 

almost unsaturated. The population results at 0 /2it =0.1 are: 

P 



”u /n 0H 

(l-r) 

r /a 

a 

82 p 

3.0 x io~ 2 

0.1 

160 

115 P 

9.0 X 10" 3 

0 .l 6 

73 

134 ’ p 

3.2 X 10“ 3 

0.20 

62 


115 






and at 


Q /2it = O. 5 O 
P 


the excitation results in: 


X 

VV 

(1-r) 

r /a 

a 

82 

3.1 x 10~ 2 

0.3 

800 

96 p 

_2 

1.5 X 10 

0.1 

300 

115 p 

9.1 X 10~ 3 

0.2 

250 

134 P 

_3 

5.0 X 10 J 

0.1 

13^0 


For the OH emission of the Orion region this means an intensity pattern 

given by the following at fi /2it = O. 5 O: 

s 


X 

* 

C D 

• State 

( T b) 

2 

I^(w/cm Hz str) 

82 p 

32 

P.S. a/R = .9 

6 

1.3 x 10 

-13 

5.3 X 10 XJ 

96 p 

21 

P.S. a/R-= 0.9 

5.7 x io 5 

-13 

1.7 X 10 J 

115 H 

12 

j P.S. a/R = 0.1 

' I. 

2.6 x 10 

5.4 x 10- 15 

i '34 p 

6.2 

F.S. K = 1 

3 

1.0 x 10 

1.5 x 10“ 


where we have used the values of: 


r i4 

physical size = 6 X 10 cm 
n 0H = 50 Cm " 3 


For the w3 OH source the intensity pattern is: 
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-13 

2.9 x 10 0 
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4 o 

P.S. a/R = .6 

1.5 x 10^ 

'4.5 X 10" 13 

115 p 

24- 

F.S. K = 2 

2.8 X 10 5 

-l4 

5.9 X 10 

134 p 

20 

P.S. a/R = 0.1 

_ 4 

1.3 x 10 

-15 

2.0 X 10 ^ 


ll6 



physical size = 6 X 10 J cm 
-3 

n = • 10. cm 
On 

2 -3 

In these last two calculations, values of- n^ below 10 cm have 
been used to show that detectable IR emission can be produced even for 
low OH densities if the pumping is s’trong enough. 

D. The Detectability Limit for OH Emission 

The formula of Chapter IV can be used to determine the emission 
intensities of the OH sources in engineering units, or the characteristics 

i 

of potential infrared instruments can be used to estimate the detecta- 
bility limit in terms of the average brightness temperature of OH sources. 
Most of the current 'infrared instruments are designed as wideband detec- 
tors, and are mediocre tools for detecting narrow line emission in the 
infrared spectrum. A reasonable sensitivity limit for the narrow emission 
lines of the OH masers for the one meter infrared telescope of the NASA 

Ames Research Center when used as a broadband receiver is at an effective 

*5 o 

brightness temperature of 10^ K. At this level several of the emission 
lines given in the previous section would be easily detected, but dis- 
tinguishing the line component may be difficult. The regions where 
substantial OH emission is possible are already known to have significant 
amounts of near and far infrared emission from extended dust clouds. 

Since the angular resolution of a one meter IR telescope is 20 arc seconds, 
or an order of magnitude larger than the OH sources, a significant amount 
of diffuse background radiation will also be measured. Unless care is 
taken to switch the beam between areas of similar background, the OH 
line emission will be indistinguishable from the background. The use of 
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high resolution, infrared spectrometers can ease this difficulty and also 
lower the limit of detectability. 

The detectability of an infrared source of limited size is set by 
the available integration time during the observation and the power 
collected by the instrument. The total power on the detector from a 
source smaller than the beam size with an effective brightness tempera- 
ture of T extending over a disc of D" arc seconds in diameter in an 
B 

emission line Av wide is: 


]p= 0.96 x io“ 13 ^p^ 



watts 


( VI-D-l) 


where A is the aperature of the instrument in meters. For a sensiti- 
vity of 3 X 10 watts (Hz) an integration time of 100 seconds 

will give a signal to noise of 30 if the factors in the above relation 
are all unity, using a detection band between 5 O and 300 microns. Longer 
integration times will improve this figure, but guidance stability can 

become a problem. At this figure an effective brightness temperature of 
*5 o 

Hr K can be considered a detectable signal if the emission can be 
distinguished from the background. The extended W3 OH model would have 
an angular extent of about 2 arc seconds, therefore detectable IR emission 
is expected for most pumping models. The Orion OH source has a physical 
extent that puts D" at' less than 0.1 arc seconds, thus the IR lines would 
be detectable only if the pumping was very strong, 

E. Summary and Conclusions 

A number of infrared emission lines from the OH molecule have been 
shown to undergo maser amplification to detectable levels from the Orion OH 
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and the W3 OH radio sources. The energy structure of the transition 
matrix of the OH molecule was carefully studied to determine the available 
pathways of excitation .in molecular OH. A transfer equation which was 
coupled to the complete population solution was solved for a spherical 
geometry of general saturation state for the characteristic degree of 
saturation and intensity as controlled by the physical parameters of the 
cloud and the strength of the pumping. • The transfer matrix of the OH 
molecule was then combined with the effects of pumping by ultraviolet, 
infrared^ and collisional mechanisms to determine the population distribu- 
tions and to determine the degree of population inversion occurring. 

These population solutions combined with- the transfer equation to estimate 

infrared line intensities at wavelengths near one hundred microns, and 

r To 

to predict effective brightness temperatures between 10 and 10 1 K for 
the detectable infrared lines. 

Because the microwave masers are fully saturated the effects of the 
background radiation are minimized. Since the. gain in the saturated 
cloud is a linear function of the diameter, adding a sizeable background 
brightness just allows the required size of the OH cloud to be reduced 
somewhat, or the pumping rate lowered. If amplified background radiation 
contributes one third of the emitted intensity and the background of the 
OH region in W3 or W^9 is the infrared radiation observed nearby, then the 
time correlation suggested between infi'ared and microwave emissions can 
be explained as the result of varying the background component of the 
amplified emission. Then the variations expected in the infrared emission 
will be very large. The factor for the' background intensity in the maser 
emission is nearly three orders of magnitude smaller for the infrared 
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lines than" for the microwave lines, but the infrared emission is nearly 
unsaturated, so that these small variations will be greatly amplified. 

The infrared intensities for these OH sources would thus vary over a 
range at least as large as the microwave emission, to an order of magni- 
tude greater variation. 

For all the pumping models considered for the OH emission there has 
been at least one proposed manner in which the microwave inversions 
could be produced. For pumping by radiational processes a complex propa- 
gation effect is necessary in order to produce the parity selection which 
populates one' of the lambda doublet levels over the other and hence 
creates the microwave emission. For infrared pumping the effects of 
radiation trapping in the 82 micron line are needed, and it is difficult 

l4 

to understand the needed amounts of trapping in regions only 10 to 
15 

10 cm across. The population inversion produced by a collisional pumping 
model which greatly simplifies the expected differential collision effects 
is in its simplicity a marked contrast. The collisions with moderate 
density molecules easily produced the needed microwave inversion even 
at densities which precluded detectable infrared emission, so that it is 
difficult to see how collisions cannot play an important role in the 
excitation of the OH masers. 

In some of the OH sources similar to the Orion OH source the presence 
of large amounts of dust and other obscuration could significantly 
attenuate the infrared line emission; however this process itself results 
in significant amounts of heating for the obscuring dust, and could be a 
significant factor in the processes occurring in the region such as the 
formation of condensations and protostars. The mechanism of transferring 
energy from a very hot source such as an O star to the dust is easily 
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accomplished by such line radiation. If the OH emission is pumped by 

a collisional process, on the other hand, the extraction of the kinetic 

energy from the cloud by the infrared radiation could lead the gas and 

dust evolution in a different direction. It seems then that in regions 

where significant infrared OH emission occurs that emission can be an 

important element in the local evolution. 

Figure 10a shows that for values of C D under ten there is little 

more than an order of magnitude between the unsaturated limit intensity 

and the inversion limit intensity. Since the inversion limit is the 

point at which the population difference becomes zero, it is expected 

* 

that the intensity of the emission for small enough values of C D will 
be nearly independent of the population of the lower state. This of 
course is the case of integrated spontaneous , emission, where 

, \ ( li u /n 0H> n 0H A uI hV < D > 

vv/ . Av 

with D the physical diameter. Since most of the IR emission occurs at 
* 

small values of C D this allows us a check of the maser calculated inten- 
sities. The 135 micron line of the collisional. pumping model with 

k -3 ^ 

n__ = 10 cm would have a spontaneous emission of 

H 2 

1.9 X 10 ^ w cm Hz str , or less than two orders of magnitude below 

the maser emission, as expected. That the maser amplification .only 
provides a gain of a hundred is not unreasonable for the relatively low 
populations found in the IR excited states. Further, the pumping rates, 
while high enough to prevent complete saturation, have the IR masers 
rather close to the critical core size. The exponential gain of the ' 
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maser comes only from this unsaturated core and little intensity is 
contributed by the saturated shell. Finally, this level of spontaneous 
emission is available only for the pumping described, since normal 
thermal mechanisms could not provide sufficient population to achieve 
even the above intensity. 
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